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Ever since its discovery, the Berry phase has permeated through all branches of physics. Over the 
last three decades, it was gradually realized that the Berry phase of the electronic wave function 
can have a profound effect on material properties and is responsible for a spectrum of phenomena, 
such as ferroelectricity, orbital magnetism, various (quantum/anomalous/spin) Hall effects, and 
quantum charge pumping. This progress is summarized in a pedagogical manner in this review. 
We start with a brief summary of necessary background, followed by a detailed discussion of the 
Berry phase effect in a variety of solid state applications. A common thread of the review is the 
semiclassical formulation of electron dynamics, which is a versatile tool in the study of electron 
dynamics in the presence of electromagnetic fields and more general perturbations. Finally, we 
demonstrate a re-quantization method that converts a semiclassical theory to an effective quantum 
theory. It is clear that the Berry phase should be added as a basic ingredient to our understanding 
of basic material properties. 
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I. INTRODUCTION 
A. Topical overview 

In 1984, Michael Berry wrote a paper that has gener- 
ated immense interests throughout the different fields of 
physics including quantum chemistry (Berry, 1984). This 
is about adiabatic evolution of an eigenenergy state when 
the external parameters of a quantum system change 
slowly and make up a loop in the parameter space. In 
the absence of degeneracy, the eigenstate will surely come 
back to itself when finishing the loop, but there will be a 
phase difference equal to the time integral of the energy 
(divided by h) plus an extra which is later commonly 
called the Berry phase. 

The Berry phase has three key properties that make 
the concept important. First, it is gauge invariant. The 
eigen-wavefunction is defined by a homogeneous linear 
equation (the eigenvalue equation), so one has the gauge 
freedom of multiplying it with an overall phase factor 
which can be parameter dependent. The Berry phase 
is unchanged (up to unessential integer multiple of 2tt) 
by such a phase factor, provided the eigen-wavefunction 
is kept to be single valued over the loop. This property 
makes the Berry phase physical, and the early experimen- 
tal studies were focused on measuring it directly through 
interference phenomena. 

Second, the Berry phase is geometrical. It can be writ- 
ten as a line-integral over the loop in the parameter space, 
and does not depend on the exact rate of change along 
the loop. This property makes it possible to express the 
Berry phase in terms of local geometrical quantities in 
the parameter space. Indeed, Berry himself showed that 
one can write the Berry phase as an integral of a field, 
which we now call as the Berry curvature, over a surface 
suspending the loop. A large class of applications of the 
Berry phase concept occur when the parameters them- 
selves are actually dynamical variables of slow degrees of 
freedom. The Berry curvature plays an essential role in 
the effective dynamics of these slow variables. The vast 
majority of applications considered in this review are of 
this nature. 

Third, the Berry phase has close analogies to gauge 
field theories and differential geometry (Simon, 1983). 
This makes the Berry phase a beautiful, intuitive and 
powerful unifying concept, especially valuable in today's 
ever specializing physical science. In primitive terms, 
the Berry phase is like the Aharonov-Bohm phase of a 
charge particle traversing a loop including a magnetic 
flux, while the Berry curvature is like the magnetic field. 
The integral of the Berry curvature over closed surfaces, 



such as a sphere or torus, is known to be topological and 
quantized as integers (Chern numbers). This is analo- 
gous to the Dirac monopoles of magnetic charges that 
must be quantized in order to have a consistent theory of 
quantum mechanical theory for charged particles in mag- 
netic fields. Interestingly, while the magnetic monopoles 
are yet to be detected in the real world, the topological 
Chern numbers have already found correspondence with 
the quantized Hall conductance plateaus in the spectac- 
ular quantum Hall phenomenon (Thoulcss et al, 1982). 

This review is about applications of the Berry phase 
concept in solid state physics. In this field, we are typ- 
ically interested in macroscopic phenomena which are 
slow in time and smooth in space in comparison with 
the atomic scales. Not surprisingly, the vast majority of 
applications of the Berry phase concept are found here. 
This field of physics is also extremely diverse, and we 
have many layers of theoretical frameworks with differ- 
ent degrees of transparency and validity (Bohm et al, 
2003; Shapcrc and Wilczek, 1989a). Therefore, a unify- 
ing organizing principle such as the Berry phase concept 
is particularly valuable. 

We will focus our attention on electronic properties, 
which play a dominant role in various aspects of mate- 
rial properties. The electrons are the glue of materials 
and they are also the agents responding swiftly to exter- 
nal fields and giving rise to strong and useful signals. A 
basic paradigm of theoretical framework is based on the 
assumption that electrons are in Bloch waves traveling 
more or less independently in periodic potentials of the 
lattice, except that the Pauli exclusion principle has to be 
satisfied and electron-electron interactions are taken care 
of in some self-consistent manner. Much of our intuition 
on electron transport are derived from the semiclassical 
picture that the electrons behave almost as free particles 
in response to external fields provided one uses the band 
energy in place of the free-particle dispersion. Partly for 
this reason, first-principles studies of electronic proper- 
ties tend to only document the energy band structures 
and various density profiles. 

There have been overwhelming evidences that such 
a simple picture cannot give complete account of ef- 
fects to first order in the fields. The anomalous ve- 
locity is a prime example, whose existence were pre- 
dicted in the early 50s (Karplus and Luttinger, 1954) and 
whose experimental relevance were established only re- 
cently (Lee et al, 2004a, b). The usual derivation of the 
anomalous velocity is based on a linear response analy- 
sis: the velocity operator has off-diagonal elements and 
electric field mixes the bands, so that the expectation 
value of the velocity acquires an additional term propor- 
tional to the field other than the usual group velocity 
in the original band (Adams and Blount, 1959; Blount, 
1962b). The anomalous velocity were first recognized in 
mid and late 90s as an effect of the Berry curvature intrin- 
sic to each of the bands, whose existence has really noth- 
ing to do with the external field (Chang and Niu, 1995, 
1996; Sundaram and Niu, 1999). This understanding en- 
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abled us to make a direct connection with the topologi- 
cal Chern number formulation of the quantum Hall effect 
(Kohmoto, 1985; Thoulcss et ai, 1982), giving incentive 
as well as confidence in our pursuit of the eventually suc- 
cessful intrinsic explanation of the anomalous Hall effect 
(Fang et ai, 2003; Jungwirth et ai, 2002; Taguchi et ai, 
2001; Yao et ai, 2004; Zeng et ai, 2006). . 

Interestingly enough, the traditional view cannot even 
explain some basic effects to zeroth order of the fields. 
The two basic electromagnetic properties of solids as a 
medium are the electric polarization and magnetization, 
which can exist in the absence of electric and magnetic 
fields in ferroelectric and ferromagnetic materials. Their 
classic definition were based on the picture of bound 
charges and currents, but these are clearly inadequate 
for the electronic contribution and it was known that 
the polarization and orbital magnetization cannot be de- 
termined from the charge and current densities in the 
bulk of a crystal at all. A breakthrough on electric 
polarization were made in early 90s by linking it with 
the phenomenon of adiabatic charge transport and ex- 
pressing it in terms of the Berry phase 1 across the 
entire Brillouin zone (King-Smith and Vanderbilt, 1993; 
Resta, 1992; Vanderbilt and King-Smith, 1993). Based 
on the Berry phase formula, one can now routinely calcu- 
late polarization related properties using first principles 
methods, with a typical precision of the density func- 
tional theory. The breakthrough on orbital magnetiza- 
tion came only recently, showing that it not only consists 
of the orbital moments of the quasi particles but also con- 
tains a Berry curvature contribution of topological ori- 
gin (Shi et ai, 2007; Thonhauscr et al., 2005; Xiao et ai, 
2005). 

In this article, we will follow the traditional semiclas- 
sical formalism of quasiparticle dynamics, only to make 
it more rigorous by including the Berry curvatures in 
the various facets of the phase space including the adia- 
batic time parameter. All of the above mentioned effects 
are transparently revealed with complete precision of the 
full quantum theory. Related effects on anomalous ther- 
moelectric, valley Hall and magneto-transport are easily 
predicted, and other effects due to crystal deformation 
and order parameter inhomogeneity can also be explored 
without difficulty. Moreover, by including a side-jump 
term (which is itself a kind of Berry phase effect) into 
the usual Boltzmann equation and associated transport 
during collisions as well as anomalous transport between 
collisions, the scmiclassical theory can also reproduce all 
the intricacies of linear response theory in weakly disor- 
dered systems (Sinitsyn, 2008). On a microscopic level, 
although the electron wavepacket dynamics is yet to be 
directly observed in solids, the formalism has been repli- 



Also called Zak's phase, it is independent of the Berry curva- 
ture which only characterize Berry phases over loops continu- 
ously shrinkable to zero (Zak, 1989a). 



cated for light transport in photonic crystals, where the 
associated Berry phase effects are vividly exhibited in 
experiments (Bliokh et ai, 2008). Finally, it is possi- 
ble to generalize the semiclassical dynamics in a single 
band into one with degenerate or nearly degenerate bands 
(Culccr et ai, 2005; Shindou and Imura, 2005), and one 
can study transport phenomena where interband coher- 
ence effects such as in spin transport, only to realize 
that the Berry curvatures and quasiparticle magnetic mo- 
ments become non-abealian (i.e., matrices). 

The semiclassical formalism is certainly amendable to 
include quantization effects. For integrable dynamics, 
such as Bloch oscillations and cyclotron orbits, one can 
use the Bohr-Sommerfeld or EKB quantization rule. The 
Berry phase enters naturally as a shift to the classi- 
cal action, affecting the energies of the quantized levels, 
e.g., the Wannier-Stark ladders and the Landau levels. 
A high point of this kind of applications is the expla- 
nation of the intricate fractal-like Hofstadter spectrum 
(Chang and Niu, 1995, 1996). A recent breakthrough 
has also enabled us to find the density of quantum states 
in the phase space for the general case (including non- 
integrable systems) (Xiao et ai, 2005), revealing Berry- 
curvature corrections which should have broad impacts 
on calculations of equilibrium as well as transport prop- 
erties. Finally, one can execute a generalized Peierls sub- 
stitution relating the physical variables to the canon- 
ical variables, turning the semiclassical dynamics into 
a full quantum theory valid to first order in the fields 
(Chang and Niu, 2008). Spin-orbit coupling and various 
mysterious Yafet terms are all found a simple explanation 
from this generalized Peierls substitution. 

Therefore, it is clear that Berry phase effects in solid 
state physics are not something just nice to be found here 
and there, the concept is essential for a coherent under- 
standing of all the basic phenomena. It is the purpose of 
this review to summarize a theoretical framework which 
continues the traditional semiclassical point of view but 
with a much broader range of validity. It is necessary 
and sufficient to include the Berry phases and gradient 
energy corrections in addition to the energy dispersions 
in order to account all phenomena up to first order in the 
fields. 



B. Organization of the review 

This review can be divided into three main parts. In 
sec. II we consider the simplest example of Berry phase in 
crystals: the adiabatic transport in a band insulator. In 
particular, we show that induced adiabatic current due 
to a time-dependent perturbation can be written as a 
Berry phase of the electronic wave functions. Based on 
this understanding, the modern theory of electric polar- 
ization is reviewed. In sec. Ill the electron dynamics in 
the presence of an electric field is discussed as a specific 
example of the time-dependent problem, and the basic 
formula developed in Sec. II can be directly applied. In 
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this case, the Berry phase manifest as transverse velocity 
of the electrons, which gives rise to a Hall current. We 
then apply this formula to study the quantum, anoma- 
lous, and valley Hall effect. 

To study the electron dynamics under spatial- 
dependent perturbations, we turn to the semiclassical 
formalism of Bloch electron dynamics, which has proven 
to be a powerful tool to investigate the influence of slowly 
varying perturbations on the electron dynamics. Sec. IV, 
we discuss the construction of the electron wave packet 
and show that the wave packet carries an orbital mag- 
netic moment. Two applications of the wave packet ap- 
proach, the orbital magnetization, and anomalous ther- 
moelectric transport in ferromagnet are discussed. In 
Sec. V the wave packet dynamics in the presence of elec- 
tromagnetic fields is studied. We show that the Berry 
phase not only affects the equations of motion, but also 
modifies the electron density of states in the phase space, 
which can be changed by applying a magnetic field. The 
formula of orbital magnetization is rederived using the 
modified density of states. We also presented a compre- 
hensive study of the magneto-transport in the presence 
of the Berry phase. The electron dynamics under more 
general perturbations is discussed in Sec. VI. We again 
present two applications: electron dynamics in deformed 
crystals and polarization induced by inhomogeneity. 

In the remaining part of the review, we focus on the re- 
quantization of the semiclassical formulation. In Sec. VII, 
the Bohr-Sommcrfcld quantization is reviewed in details. 
With its help, one can incorporate the Berry phase effect 
into the Wannier-Stark ladders and the Landau levels 
very easily. In Sec. VIII, we show that the same semi- 
classical approach can be applied to systems subject to a 
very strong magnetic field. One only has to separate the 
field into a quantization part and a perturbation. The 
former should be treated quantum mechanically to ob- 
tain the magnetic Bloch band spectrum while the latter 
is treated perturbatively. Using this formalism, we dis- 
cuss the cyclotron motion, the splitting into magnetic 
subbands, and the quantum Hall effect. In Sec. IX, we 
review recent advance on the non-Abelian Berry phase 
in degenerate bands. We show that the Berry connec- 
tion now plays an explicit role in spin dynamics and is 
deeply related to the spin-orbit interaction. We then cite 
the relativistic Dirac electrons and the Kane model in 
semiconductors as two examples of application. Finally, 
we briefly discuss the re-quantization of the semiclassical 
theory and the hierarchy of effective quantum theories. 

We do not attempt to cover all of the Berry phase ef- 
fects in this review. Interested readers can consult the 
following books or review articles for many more left un- 
mentioned: Bohm et al. (2003); Chang and Niu (2008); 
Nenciu (1991); Resta (1994, 2000); Shapere and Wilczek 
(1989b); Tcufcl (2003); Thouless (1998). In this review, 
we focus on a pedagogical and self-contained approach, 
with the main machinery being the semiclassical formal- 
ism of Bloch electron dynamics (Chang and Niu, 1995, 
1996; Sundaram and Niu, 1999). We shall start with 



the simplest case, the gradually expand the formalism 
as more complicated physical situations are considered. 
Whenever a new ingredient is added, a few applications 
is provided to demonstrate the basic ideas. The vast 
number of application we discussed is a reflection of the 
universality of the Berry phase effect. 

C. Basic Concepts of The Berry phase 

In this subsection we introduce the basic concepts of 
the Berry phase. Following Berry's original paper (Berry, 
1984), we first discuss how the Berry phase arises as a 
generic feature of the adiabatic evolution of a quantum 
state. We then introduce the local description of the 
Berry phase in terms of the Berry curvature. A two- 
level model is used to demonstrate these concepts as well 
as some important properties, such as the quantization 
of the Berry phase. Our aim is to provide a minimal 
but self-contained introduction. For a detailed account 
of the Berry phase, including its mathematical founda- 
tion and applications in a wide range of fields in physics, 
we refer the readers to the books by Bohm et al. (2003); 
Shapere and Wilczek (1989b) and references therein. 

1. Cyclic adiabatic evolution 

Let us consider a physical system described by a Hamil- 
tonian that depends on time through a set of parameters, 
denoted by R = (i?i, R 2 , . . . ), i.e., 

H = H(R), R=R(t). (1.1) 

We are interested in the adiabatic evolution of the system 
as R(t) moves slowly along a path C in the parameter 
space. For this purpose, it will be useful to introduce an 
instantaneous orthonormal basis from the eigenstates of 
H(R) at each value of the parameter R, i.e., 

H(R)\n(R)) = e n (R)\n(R)) . (1.2) 

However, Eq. (1.2) alone does not completely determine 
the basis function \n(R)); it still allows an arbitrary R- 
dependent phase factor of \n(R)). One can make a phase 
choice, also known as a gauge, to remove this arbitrari- 
ness. Here we require that the phase of the basis function 
is smooth and single- valued along the path C in the pa- 
rameter space. 2 

According to the quantum adiabatic theorem (Kato, 
1950; Messiah, 1962), a system initially in one of its 



2 Strictly speaking, such a phase choice is guaranteed only in fi- 
nite neighborhoods of the parameter space. In the general case, 
one can proceed by dividing the path into several such neigh- 
borhoods overlapping with each other, then use the fact that in 
the overlapping region the wave functions arc related by a gauge 
transformation of the form (1.7). 
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eigenstates \n(R(0))) will stay as an instantaneous eigen- 
state of the Hamiltonian H(R(t)) throughout the pro- 
cess. Therefore the only degree of freedom we have is the 
phase of the quantum state. Write the state at time t as 

\^ n (t)) = e^We-^o'^^ln^t))) , (1.3) 

where the second exponential is known as the dynamical 
phase factor. Inserting Eq. (1.3) into the time-dependent 
Schrodinger equation 



ih^\Mt)) = H(R(t))\Mt)) 



(1.4) 



and multiplying it from the left by (n(R(t))\, one finds 
that 7„ can be expressed as a path integral in the param- 
eter space 



7n = J^dR-An{R) 



where A^i(R) is a vector- valued function 
A n (R) = i(n(R)\~\n(R)) 



(1.5) 



(1.6) 



This vector A^ l (R) is called the Berry connection or the 
Berry vector potential. Equation (1.5) shows that in ad- 
dition to the dynamical phase, the quantum state will 
acquire an additional phase 7„ during the adiabatic evo- 
lution. 

Obviously, An(R) is gauge-dependent. If we make a 
gauge transformation 



i(R))^e iaR) \n(R)) 



(1.7) 



This shows that 7„ can be only changed by an integer 
multiple of 2n under the gauge transformation (1.7) and 
it cannot be removed. Therefore for a closed path, j n 
becomes a gauge-invariant physical quantity, now known 
as the Berry phase or geometric phase in general; it is 
given by 
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= j> dR ■ A n {R) ■ 



(1.10) 



From the above definition, we can see that the Berry 
phase only depends on the geometric aspect of the closed 
path, and is independent of how R(t) varies in time. The 
explicit time-dependence is thus not essential in the de- 
scription of the Berry phase and will be dropped in the 
following discussion. 



2. Berry curvature 

It is useful to define, in analogy to electrodynamics, a 
gauge field tensor derived from the Berry vector poten- 
tial: 



%u{R) = ^A n u {R)- 1 ^A;iRA 



dR" uy ' dR 
dn(R) dn(R) 



dR» 1 dR" 



) - {v <-> n) 



(1.11) 



This field is called the Berry curvature. Then according 
to Stokes's theorem the Berry phase can be written as a 
surface integral 



1 



7„ = / dR" A dR v -n^(R) , 



(1.12) 



with C(-R) being an arbitrary smooth function, A*n(R) 
transforms according to 

An(R) - An(R) - ^C(R) • (1-8) 

Consequently, the phase 7„ given by Eq. (1.5) will be 
changed by £(K(0)) — ((R(T)) after the transformation, 
where H(0) and R(T) are the initial and final points 
of the path C. This observation has led Fock (1928) to 
conclude that one can always choose a suitable ((R) such 
that 7„ accumulated along the path C is canceled out, 
leaving Eq. (1.3) with only the dynamical phase. Because 
of this, the phase j n has long been deemed unimportant 
and it was usually neglected in the theoretical treatment 
of time-dependent problems. 

This conclusion remained unchallenged until Berry 
(1984) reconsidered the cyclic evolution of the system 
along a closed path C with R(T) = R(0). The phase 
choice we made earlier on the basis function \n(R)) re- 
quires in the gauge transformation, Eq. (1.7), to 
be single-valued, which implies 



where S is an arbitrary surface enclosed by the path C. 
It can be verified from Eq. (1.11) that, unlike the Berry 
vector potential, the Berry curvature is gauge invariant 
and thus observable. 

If the parameter space is three-dimensional, Eqs. (1.11) 
and (1.12) can be recasted into a vector form 

n n (R) = V R X An{R) , (1.11') 

7 „= / dS-n n (R) . (1.12') 
Js 

The Berry curvature tensor fi™„ and vector r2„ is related 
by fi™ = e MI ,{(r2„)£ with being the Levi-Civita an- 
tisymmetry tensor. The vector form gives us an intuitive 
picture of the Berry curvature: it can be viewed as the 
magnetic field in the parameter space. 

Besides the differential formula given in Eq. (1.11), the 
Berry curvature can be also written as a summation over 
the eigenstates: 



{n\§^\n'){n'\^\n)-{v^n) 



C(R(0)) - C(R(T)) = 2tt x integer 



(1.9) 



(e r , 



(1.13) 
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It can be obtained from Eq. (1.11) by using the relation 
(n\dH/dR\n') = (dn/dR\n')(e n - e n >) for ri ^ n. The 
summation formula has the advantage that no differenti- 
ation on the wave function is involved, therefore it can be 
evaluated under any gauge choice. This property is par- 
ticularly useful for numerical calculations, in which the 
condition of a smooth phase choice of the eigenstates is 
not guaranteed in standard diagonalization algorithms. 
It has been used to evaluate the Berry curvature in crys- 
tals with the eigenfunctions supplied from first-principles 
calculations (Fang et al., 2003; Yao et al., 2004). 

Equation (1.13) offers further insight on the origin of 
the Berry curvature. The adiabatic approximation we 
adopted earlier is essentially a projection operation, i.e., 
the dynamics of the system is restricted to the nth energy 
level. In view of Eq. (1.13), the Berry curvature can be 
regarded as the result of the "residual" interaction of 
those projected-out energy levels. In fact, if all energy 
levels are included, it follows from Eq. (1.13) that the 
total Berry curvature vanishes for each value of R, 

Y / ^u( R ) = °- ( L14 ) 

n 

This is the local conservation law for the Berry curva- 
ture. 3 Equation (1.13) also shows that fi"„(.R) be- 
comes singular if two energy levels e n (R) and e n '(R) 
are brought together at certain value of R. This degen- 
eracy point corresponds to a monopolc in the parameter 
space; an explicit example is given below. If the degen- 
erate points form a string in the parameter space, it is 
known as the Dirac string. 

So far we have discussed the situation where a sin- 
gle energy level can be separated out in the adia- 
batic evolution. However, if the energy levels are 
degenerate, then the dynamics must be projected 
to a subspace spanned by those degenerate eigen- 
states. Wilczek and Zee (1984) showed that in this 
situation non-Abelian Berry curvature naturally arises. 
Culccr et al. (2005); Shindou and Imura (2005) have dis- 
cussed the non-Abelian Berry curvature in the context of 
degenerate Bloch bands. In the following we shall focus 
on the Abelian formulation and defer the discussion of 
the non-Abelian Berry curvature to Sec. IX. 

Compared to the Berry phase which is always asso- 
ciated with a closed path, the Berry curvature is truly 
a local quantity. It provides a local description of the 
geometric properties of the parameter space. More 
importantly, just like a magnetic field can affect the 
electron dynamics, the Berry curvature also directly 



3 The conservation law is obtained under the condition that the 
full Hamiltonian is known. However, in practice one usually deals 
with effective Hamiltonians which are obtained through some 
projection process of the full Hamiltonian. Therefore there will 
always be some "residual" Berry curvature accompanying the ef- 
fective Hamiltonian. See Chang and Niu (2008) and discussions 
in Sec. IX. 



participates in the dynamics of the adiabatic parame- 
ters (Kuratsuji and Iida, 1985). In this sense, the Berry 
curvature is a more fundamental quantity than the Berry 
phase. 



3. Example: The two-level system 

Let us consider a concrete example: a two-level system. 
The purpose to study this system is two- fold. Firstly, as 
a simple model, it demonstrates the basic concepts as 
well as several important properties of the Berry phase. 
Secondly, it will be repeatedly used later in this article, 
although in different physical context. It is therefore use- 
ful to go through the basis of this model. 

The generic Hamiltonian of a two-level system takes 
the following form 

H = h(R)-a, (1.15) 

where <x is the Pauli matrices. Despite its simple 
form, the above Hamiltonian describes a number of 
physical systems in condensed matter physics for which 
the Berry phase effect has been discussed. Exam- 
ples include spin-orbit coupled systems (Culcer et al, 
2003; Liu et al., 2008), linearly conjugated diatomic 
polymers (Rice and Mele, 1982; Su et al., 1979), 
one-dimensional ferroelectrics (Onoda et al, 2004b; 
Vanderbilt and King-Smith, 1993), graphene (Haldanc, 
1988; Semenoff, 1984), and Bogoliubov quasiparti- 
cles (Zhang et al., 2006). 

Parameterize h by its azimuthal angle 9 and polar an- 
gle (j>, h = h(sin 9 cos 4>, sin 6 sin <fi, cos 9). The two eigen- 
states with energies ±/i is 

s /sinfe"^\ , s /cos|e-^\ , 

«- = 2 e . K> = e ■ !-16 

V - cos I / V sin | / 

We are, of course, free to add an arbitrary phase to these 
wave functions. Let us consider the lower energy level. 
The Berry connection is given by 

A e = (u\id e u) = , (1.17a) 

Q 

Aj, = {u\id^u) = sin 2 - , (1.17b) 

and the Berry curvature is 

Q. e <t> = deA^ - d^A g = ^sin0 . (1.18) 

However, the phase of \u-) is not defined at the south 
pole (6 = 7r). We can choose another gauge by multiply- 
ing by e 1 ^ so that the wave function is smooth and 
single valued everywhere except at the north pole. Under 
this gauge we find Ae = and A<j> = — cos 2 § , and the 
same expression for the Berry curvature as in Eq. (1.18). 
This is not surprising because the Berry curvature is a 
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gauge-independent quantity and the Berry connection is 
not. 4 

If h(R) depends on a set of parameters R, then 
1 d((f>, cos 9) 



2d(R 1 ,R 2 ) ■ 



(1.19) 



Several important properties of the Berry curvature 
can be revealed by considering the specific case of h = 
(x,y,z). Using Eq. (1.19), we find the Berry curvature 
in its vector form 



~ 1 h 

ft = 

2/i 3 



(1.20) 



One recognizes that Eq. (1.20) is the field generated by 
a monopole at the origin h = (Dirac, 1931; Sakurai, 
1993; Wu and Yang, 1975), where the two energy levels 
become degenerate. Therefore the degeneracy points act 
as sources and drains of the Berry curvature flux. Inte- 
grate the Berry curvature over a sphere containing the 
monopole, which is the Berry phase on the sphere; we 
find 



1 

2ir j S 2 



£led> = 1 



(1.21) 



In general, the Berry curvature integrated over a closed 
manifold is quantized in the units of 2tt and equals to the 
net number of monopoles inside. This number is called 
the Chern number and is responsible for a number of 
quantization effects discussed below. 



D. Berry phase in Bloch bands 

In the above we have introduced the basic concepts 
of the Berry phase for a generic system described by a 
parameter-dependent Hamiltonian. We now consider its 
realization in crystalline solids. As we shall see, the band 
structure of crystals provides a natural platform to inves- 
tigate the occurrence of the Berry phase effect. 

Within the independent electron approximation, the 
band structure of a crystal is determined by the following 
Hamiltonian for a single electron: 



H 



V 

2m 



V{r) 



(1.22) 



where V(r + a) = V(r) is the periodic potential with 
a being the Bravais lattice vector. According to Bloch's 
theorem, the eigenstates of a periodic Hamiltonian satisfy 
the following boundary condition 



ip nq (r + a) 



(1.23) 



4 Another way to calculate the Berry curvature is to use Eq. (1.13) 
directly. 



where n is the band index and hq is the crystal momen- 
tum, which resides in the Brillouin zone. Thus the sys- 
tem is described by a q-independent Hamiltonian with 
a (/-dependent boundary condition, Eq. (1.23). To com- 
ply with the general formalism of the Berry phase, we 
make the following unitary transformation to obtain a 
q-dependent Hamiltonian: 



H(q) 



2m 



V(r) . (1.24) 



The transformed eigenstate u nq (r) — e^ lq ' v ip nq (r) is just 
the cell-periodic part of the Bloch function. It satisfies 
the strict periodic boundary condition 



(r- 



•'"II q 



(r) 



(1.25) 



This boundary condition ensures that all the eigenstates 
live in the same Hilbert space. We can thus identify the 
Brillouin zone as the parameter space of the transformed 
Hamiltonian H(q), and \u n (q)) as the basis function. 

Since the q-dependence of the basis function is inherent 
to the Bloch problem, various Berry phase effects are 
expected in crystals. For example, if q is forced to vary 
in the momentum space, then the Bloch state will pick 
up a Berry phase: 



7„ = fdq 



(u„(q)\iV q \u n (q)) 



(1.26) 



We emphasize that the path C must be closed to make 
7„ a gauge-invariant quantity with physical significance. 

Generally speaking, there are two ways to gener- 
ate a closed path in the momentum space. One 
can apply a magnetic field, which induces a cy- 
clotron motion along a closed orbit in the q-space. 
This way the Berry phase can manifest in various 
magneto-oscillatory effects (Mikitik and Sharlai, 1999, 
2004, 2007), which have been observed in metallic com- 
pound LaRhIn 5 (Goodrich et ai, 2002), and most re- 
cently, graphene systems (Novoselov et ai, 2005, 2006; 
Zhang et ai, 2005). Such a closed orbit is possible only in 
two or three-dimensional systems (see Sec. VILA). Fol- 
lowing our discussion in Sec. I.C, we can define the Berry 
curvature of the energy bands, given by 



« n (g) = V, 



x (U, 



,(q)\iV q \u n (q)) 



(1.27) 



The Berry curvature S~2 ?l (q) is an intrinsic property of 
the band structure because it only depends on the wave 
function. It is nonzero in a wide range of materials, in 
particular, crystals with broken time-reversal or inversion 
symmetry. In fact, once we have introduced the concept 
of the Berry curvature, a closed loop is not necessary be- 
cause the Berry curvature itself is a local gauge-invariant 
quantity. It is now well recognized that information of 
the Berry curvature is essential in a proper description of 
the dynamics of Bloch electrons, which has various effects 
on transport and thermodynamic properties of crystals. 

One can also apply an electric field to cause a linear 
variation of q. In this closed path is realized when 
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q sweeps the entire Brillouin zone. To see this, we note 
that the Brillouin zone has the topology of a torus: the 
two points q and q + G can be identified as the same 
point, where G is the reciprocal lattice vector. This can 
be seen by noting that \ip n (q)) and \ip n (q + G)) satisfy 
the same boundary condition in Eq. (1.23), therefore they 
can at most differ by a phase factor. The torus topology 
is realized by making the phase choice \ip n (q)) — \4>n(q + 
G)). Consequently, \u n (q)) and \u n (q + G)) satisfy the 
following phase relation 



.(*)> = 



JGr 



u n (q + G)) 



(1.28) 



This gauge choice is called the periodic gauge (Rcsta, 
2000). 

In this case, the Berry phase across the Brillouin zone 
is called Zak's phase (Zak, 1989a) 



In 



BZ 



dq ■ {u n (q)\iV q \u n (q)) . 



(1.29) 



It plays an important role in the formation of Wannier- 
Stark ladders (Wannier, 1962) (see Sec. VII. B). We note 
that this phase is entirely due to the torus topology of the 
Brillouin zone, and it is the only way to realize a closed 
path in a one-dimensional parameter space. By analyz- 
ing the symmetry properties of Wannier functions (Kohn, 
1959) of a one-dimensional crystal, Zak (1989a) showed 
that 7„ is either or ir in the presence of inversion 
symmetry; a simple argument is given in Sec. II. C. If 
the crystal lacks inversion symmetry, 7„ can assume any 
value. Zak's phase is also related to macroscopic polariza- 
tion of an insulator (King-Smith and Vanderbilt, 1993; 
Resta, 1994; Sipe and Zak, 1999) (see Sec. II.C). 



II. ADIABATIC TRANSPORT AND ELECTRIC 
POLARIZATION 

One of the earlier examples of the Berry phase ef- 
fect in crystals is the adiabatic transport in a one- 
dimensional band insulator, first considered by Thouless 
(1983). He found that if the potential varies slowly in 
time and returns to itself after some time, the particle 
transport during the time cycle can be expressed as a 
Berry phase and it is an integer. This idea was later 
generalized to many-body systems with interactions and 
disorder, provided that the Fermi energy always lies in 
a bulk energy gap during the cycle (Niu and Thouless, 
1984). Avron et al. (1988) studied the adiabatic trans- 
port in multiply connected systems. The scheme of adi- 
abatic transport under one or several controlling param- 
eters has proven very powerful. It opened the door 
to the field of parametric charge pumping (Brouwer, 
1998; Niu, 1990; Switkes et al, 1999; Talyanskii et al, 
1997; Zhou et al., 1999). It also provides a firm foun- 
dation to the modern theory of polarization devel- 
oped in the early 90 's (King-Smith and Vanderbilt, 1993; 
Ortiz and Martin, 1994; Resta, 1994). 



A. Adiabatic current 

Let us consider a one-dimensional band insulator un- 
der a slowly varying time-dependent perturbation. We 
assume the perturbation is periodic in time, i.e., the 
Hamiltonian satisfies 



H(t + T) = H(t) . 



(2.1) 



Since the time-dependent Hamiltonian still has the trans- 
lational symmetry of the crystal, its instantaneous eigen- 
states has the Bloch form e lqx \u n (q, t)). It is convenient 
to work with the g-space representation of the Hamilto- 
nian H{q,t) [see Eq. (1.24)] with eigenstates \u n (q,t)). 
We note that under this parametrization, the wave vec- 
tor q and time t are put on an equal footing as both are 
independent cooridinates of the parameter space. 

We are interested in the adiabatic current induced by 
the variation of external potentials. Apart from an unim- 
portant overall phase factor and up to first order in the 
rate of the change of the Hamiltonian, the wave function 
is given by 



\u n ) 



(2.2) 



The velocity operator in the g-representaion has the form 
v(q,t) = dH(q,t)/d(Tiq). Hence, the average velocity in 
a state of given q is found to first order as 



v n (q) 



dsn(q) 
hdq 



Un \^\u n ,){u n ,\^) 



da ,1 



c.c 



£77/ 



}. 



(2.3) 

where c.c. denotes the complex conjugate. Using the fact 
that (u n \dH/dq\u n >) — (e„ — e n >)(du n /dq\u n i) and the 
identity Y^ n '\ u n')( u n'\ = 1, we find 



v n (q) 



de n (g) 
hdq 



r , 8u n 8u n 

~dq~\~dt 



"~dt\~dq~ 



(2.4) 



The second term is exactly the Berry curvature f2™ t de- 
fined in the parameter space (q,t) [see Eq. (1.11)]. There- 
fore the above equation can be recasted into a compact 
form 



v n (q) 



de n (q) 
hdq 



n 



qt ■ 



(2.5) 



Upon integration over the Brillouin zone, the zeroth 
order term given by the derivative of the band energy 
vanishes, and only the first order term survives. The 
induced adiabatic current is given by 



J = 



E 



BZ 



2tt qt 



(2.6) 



where the sum is over filled bands. We have thus derived 
the remarkable result that the adiabatic current induced 
by a time-dependent perturbation in a band is equal 
to the q-integral of the Berry curvature f2™ t (Thouless, 
1983). 
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D=(Q,l) C=(l,l) 




(a) 



A=(Q,0) B=(l,0) x 
(b) 



FIG. 1 (a) A torus with its surface parameterized by (q,t). 
The control parameter R(t) runs in circle along the t direc- 
tion, (b) The equivalence of a torus: a rectangle with peri- 
odic boundary conditions: AB — DC and AD = BC. To use 
Stokes's theorem, we relax the boundary condition and allow 
the wave functions on parallel sides have different phases. 



B. Quantized adiabatic particle transport 

Next we consider the particle transport for the nth 
band over a time cycle, given by 



T 

dt 

JBZ 



2tt " 



(2.7) 



Since the Hamiltonian H(q, t) is periodic in both t and q, 
the parameter space of H(q, t) is a torus, schematically 
shown in Fig. 1(a). By definition (1.12), 27rc„ is nothing 
but the Berry phase over the torus. 

In Sec. I.C. 3, we showed that the Berry phase over a 
closed manifold, the surface of a sphere S 2 in that case, is 
quantized in the unit of 2tt. Here we prove that it is also 
true in the case of a torus. Our strategy is to evaluate 
the surface integral (2.7) using Stokes's theorem, which 
requires the surface to be simply connected. To do that, 
we cut the torus open and transform it into a rectangle, as 
shown in Fig. 1(b). The basis function along the contour 
of the rectangle is assumed to be single- valued. Introduce 
x = t/T and y — q/2n. According to Eq. (1.10), the 
Berry phase in Eq. (2.7) can be written into a contour 
integral of the Berry vector potential, i.e., 



2tt 



c 



dxA x (x,0) + / dyA y (l,y) 



+ J dxA x (x,l) + J dyA v {0,y)} 



2tt 



(2.8) 



dx[A x (x,Q) - A x {xA)} 
dy[Ay(0,y)-Ay(l,y)}} , 



where the band index n is dropped for simplicity. Let us 
consider the integration over x. By definition, A x (x, y) — 
(u(x,y)\iV x \u(x,y)). Recall that |u(x, 0)) and \u(x, 1)) 
describe physically equivalent states, therefore they can 



only differ by a phase factor, i.e., e l0xtyX ^ \u(x, 1)) 
\u(x,0)). We thus have 



dx[A x {x,Q)-A x {xA)]=e x {l)-6 x {Q) . (2.9) 



Similarly, 

dy [A y (0, y) - A y (l, y)] = 0„(1) - 9 y (0) , (2.10) 
where e 19 "^ \u(y, 1)) = \u(y,0)). The total integral is 

c = - MO) + y (o) - Mi)] ■ (2-ii) 

On the other hand, using the phase matching relations 
at the four corners A, B, C, and D, 

e ie »WK0 ) l)) = K0,Q)) ) 
e <e »«Kl ) l)) = Kl,0)) ) 
e^(°)| M (l,0)) = K0,0)), 
e^«Kl,l)) = K0,l)), 

we obtain 

K0, 0)) = eW-UM-COHMOM.d)] | u ( 0) )) . (2.12) 

The single- valuedness of \u) requires that the exponent 
must be an integer multiple of 2tt. Therefore the trans- 
ported particle number c, given in Eq. (2.11), must be 
quantized. This integer is called the first Chern number, 
which characterizes the topological structure of the map- 
ping from the parameter space (g, t) to the Bloch states 
\u(q,t)). Note that in our proof, we made no reference 
to the original physical system; the quantization of the 
Chern number is always true as long as the Hamiltonian 
is periodic in both parameters. 

An intuitive picture of the quantized particle trans- 
port is the following. If the periodic potential slides its 
position without changing its shape, we expect that the 
electrons simply follow the potential. If the potential 
shifts one spatial period in the time cycle, the particle 
transport should be equal to the number of filled Bloch 
bands (double if the spin degeneracy is counted). This 
follows from the fact that there is on average one state 
per unit cell in each filled band. 



1. Conditions for nonzero particle transport for cyclic motion 

We have shown that the adiabatic particle transport 
over a time period takes the form of the Chern number 
and it is quantized. However, the exact quantization does 
not gurantcc that the electrons will be transported at the 
end of the cycle because zero is also an integer. According 
to the discussion in Sec. I.C. 3, the Chern number counts 
the net number of monopoles enclosed by the surface. 
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Therefore the number of transported electrons can be re- 
lated to the number of monopoles, which are degeneracy 
points in the parameter space. 

To formulate this problem, we let the Hamiltonian de- 
pend on time through a set of control parameters R(t), 
i.e., 

H(q,t) = H(q,R(t)) , R(t + T) = R(t) . (2.13) 
The particle transport is now given by, in terms of R, 



with the strict periodic boundary condition 



1 

2^ 



dR a / dqfl 



BZ 



(2.14) 



If R(t) is a smooth function of t, as it is usually the 
case for physical quantities, then R must have at least 
two components, say R\ and R%. Otherwise the trajec- 
tory of R[t) cannot trace out a circle on the torus [see 
Fig. 1(a)]. To find the monopoles inside the torus, we 
now relax the constraint that R\ and R2 can only move 
on the surface and extend their domains inside the torus 
such that the parameter space of (g, Ri, R2) becomes a 
toroid. Thus, the criterion for c n to be nonzero is that a 
degeneracy point must occur somewhere inside the torus 
as one varies q, R\ and i?2- In the context of quasi one- 
dimensional ferroclcctrics, Onoda et al. (2004b) have dis- 
cussed the situation where R has three components, and 
showed how the topological structure in the R space af- 
fects the particle transport. 



2. Many-body interactions and disorder 

In the above we have only considered band insulators 
of non-interacting electrons. However, in real materi- 
als both many-body interactions and disorder are ubiqui- 
tous. Niu and Thoulcss (1984) studied this problem and 
showed that in the general case the quantization of par- 
ticle transport is still valid as long as the system remains 
an insulator during the whole process. 

Let us consider a time-dependent Hamiltonian of an 
iV-particle system 

iV ~2 N 

H ( t )=Yl[§ n + U(xi,t)]+J2V(x i -x j ), (2.15) 

i i>j 

where the one-particle potential U(xi,t) varies slowly in 
time and repeats itself in period T. Note that we have not 
assumed any specific periodicity of the potential U(xi,t). 
The trick is to use the so-called twisted boundary con- 
dition by requiring that the many-body wave function 
satisfies 



$(xi 



+ L, . . .,x N ) 



i ^ii • • • i *^iv) ) 

(2-16) 

where L is the size of the system. This is equivalent to 
solving a K-dependent Hamiltonian 

H(K,t) = exp(in Xj)H(t) exp(— Xj) (2.17) 



<&(k; Xi 



:+L, ... ,x N ) = <J> k (k;xi, 



Xi, ■ ■ ■ , Xn) ■ 

(2-18) 

The Hamiltonian H (k, t) together with the boundary 
condition (2.18) describes a one-dimensional system 
placed on a ring of length L and threaded by a magnetic 
flux of (Ti/e)nL (Kohn, 1964). We note that the above 
transformation (2.17) with the boundary condition (2.18) 
is very similar to that of the one-particle case, given by 
Eqs. (1.24) and (1.25). 

One can verify that the current operator is given by 
dH(n,t)/d(hn). For each k, we can repeat the same 
steps in Sec. II. A by identifying \u n ) in Eq. (2.2) as the 
many-body ground-state |$o) an d \u n ') as the excited 
state. We have 



i(«) 



Tidn 
Tidn 



v 8k 1 dt 1 



x dt 1 8k' 



(2.19) 



ft* 



So far the derivation is formal and we still cannot see 
why the particle transport should be quantized. The key 
step is achieved by realizing that if the Fermi energy lies 
in a gap, then the current j(n) should be insensitive to the 
boundary condition specified by k (Niu and Thouless, 
1984; Thoulcss, 1981). Consequently we can take the 
thermodynamic limit and average j(n) over different 
boundary conditions. Note that k and k + 2tt/L describe 
the same boundary condition in Eq. (2.16). Therefore 
the parameter space for k and t is a torus T 2 : {0 < k < 
2ir/L,0 < t < T}. The particle transport is given by 



1 

"27 



dt 



2ir/L 



dn fL 



(2.20) 



which, according to the previous discussion, is quantized. 

We emphasize that the quantization of the particle 
transport only depends on two conditions: 

1. The ground state is separated from the excited 
states in the bulk by a finite energy gap; 

2. The ground state is non-degenerate. 

The exact quantization of the Chern number in the pres- 
ence of many-body interactions and disorder is very re- 
markable. Usually, small perturbations to the Hamilto- 
nian results in small changes of physical quantities. How- 
ever, the fact that the Chern number must be an integer 
means that it can only be changed in a discontinous way 
and does not change at all if the perturbation is small. 
This is a general outcome of the topological invariance. 

Later we show that the same quantity also appears in 
the quantum Hall effect. The expression (2.19) of the 
induced current also provides a many-body formulation 
for adiabatic transport. 
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3. Adiabatic Pumping 

The phenomenon of adiabatic transport is sometimes 
called adiabatic pumping because it can generate a dc 
current / via periodic variations of some parameters of 
the system, i.e., 

I = ecv, (2.21) 

where c is the Chern number and v is the frequency of 
the variation. Niu (1990) suggested that the exact quan- 
tization of the adiabatic transport can be used as a stan- 
dard for charge current and proposed an experimental 
realization in nanodevices, which could serve as a charge 
pump. It was later realized in the experimental study 
of acoustoelectric current induced by a surface acoustic 
wave in a one-dimensional channel in a GaAs-Al^Gai-^ 
heterostructure (Talyanskii et al., 1997). The same idea 
has led to the proposal of a quantum spin pump in an 
antiferromagentic chain (Shindou, 2005). 

Recently, much efforts have focused on adia- 
batic pumping in mesoscopic systems (Avron et al, 
2001; Brouwer, 1998; Mucciolo et al, 2002; 
Sharma and Chamon, 2001; Zheng et al., 2003; 
Zhou et al., 1999). Experimentally, both charge 
and spin pumping have been observed in a quantum 
dot system (Switkes et al, 1999; Watson et al., 2003). 
Instead of the wave function, the central quantity in a 
mesoscopic system is the scattering matrix. Brouwer 
(1998) showed that the pumped charge over a time 
period is given by 

where m labels the contact, Xi and X2 are two exter- 
nal parameters whose trace encloses the area A in the 
parameter space, a and f3 labels the conducting chan- 
nels, and S a @ is the scattering matrix. Although the 
physical description of these open systems are dramati- 
cally different from the closed ones, the concepts of gauge 
field and geometric phase can still be applied. The in- 
tegrand in Eq. (2.22) can be thought as the Berry cur- 
vature flx ± x 2 = — 2 < 3(dx 1 u\dx 2 u) if we identify the in- 
ner product of the state vector with the matrix product. 
Zhou et al. (2003) showed the pumped charge (spin) and 
is essentially the Abclian (non-Abelian) geometric phase 
associated with scattering matrix S a p. 

C. Electric Polarization of Crystalline Solids 

Electric polarization is one of the fundamental quanti- 
ties in condensed matter physics, essential to any proper 
description of dielectric phenomena of matter. Despite 
its great importance, the theory of polarization in crys- 
tals had been plagued by the lack of a proper micro- 
scopic understanding. The main difficulty lies in the 
fact that in crystals the charge distribution is periodic in 



space, for which the electric dipole operator is not well 
defined. This difficulty is most exemplified in covalent 
solids, where the electron charges are continuously dis- 
tributed between atoms. In this case, simple integration 
over charge density would give arbitrary values depend- 
ing on the choice of the unit cell (Martin, 1972, 1974). 
It has prompted the question whether the electric polar- 
ization can be defined as a bulk property. These prob- 
lems are eventually solved by the modern theory of polar- 
ization (King-Smith and Vanderbilt, 1993; Resta, 1994), 
where it is shown that only the change in polarization 
has physical meaning and it can be quantified by using 
the Berry phase of the electronic wave function. The 
resulting Berry-phase formula has been very successful 
in first-principles studies of dielectric and ferroelectric 
materials. Resta and Vanderbilt (2007) reviewed recent 
progress in this field. 

Here we discuss the theory of polarization based on the 
concept of adiabatic transport. Their relation is revealed 
by elementary arguments from macroscopic electrostat- 
ics (Ortiz and Martin, 1994). We begin with the relation 

V • P(r) = -p(r) , (2.23) 

where P(r) is the polarization density and p(r) is the 
charge density. Coupled with the continuity equation 

^+V-i = 0, (2.24) 

Eq. (2.23) leads to 

V-(^-j) = 0. (2.25) 

Therefore up to a divergence- free part, 5 the change in 
the polarization density is given by 

AP a = [ dtj a . (2.26) 
Jo 

The above equation can be interpreted in the following 
way: The polarization difference between two states is 
given by the integrated bulk current as the system adia- 
batically evolves from the initial state at t — to the final 
state at t = T. This description implies a time-dependent 
Hamiltonian H(t), and the electric polarization can be 
regarded as "unquantized" adiabatic particle transport. 
The above interpretation is also consistent with experi- 
ments, as it is always the change of the polarization that 
has been measured (Resta and Vanderbilt, 2007). 

Obviously, the time t in the Hamiltonian can be re- 
placed by any scalar that describes the adiabatic process. 



5 The divergence-free part of the current is usually given by the 
magnetization current. In a uniform system, such current van- 
ishes identically in the bulk. Hirst (1997) gave an in-depth dis- 
cussion on the separation between polarization and magnetiza- 
tion current. 
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For example, if the process corresponds to a deformation 
of the crystal, then it makes sense to use the parame- 
ter that characterizes the atomic displacement within a 
unit cell. For general purpose, we shall assume the adi- 
abatic transformation is parameterized by a scalar X(t) 
with A(0) = and A(T) = 1. It follows from Eqs. (2.6) 
and (2.26) that 



P is either or ea/2 (modulo ea). Any other value of P 
will break the inversion symmetry. Consequently, Zak's 
phase can only take the value or tt (modulo 2ir). 

King-Smith and Vandcrbilt (1993) further showed 
that, based on Eq. (2.28), the polarization per unit cell 
can be defined as the dipole moment of the Wannier 
charge density, 



AP a = 



dX 



BZ 



d( l on 

(2ir) d 9qA 



(2.27) 



P = -e 



where d is the dimensionality of the system. 
This is the Berry-phase formula obtained by 
King-Smith and Vanderbilt (1993). 

In numerical calculations, a two-point version of 
Eq. (2.27) that only involves the integration over q is 
commonly used to reduce the computational load. It is 
obtained under the periodic gauge [see Eq. (1.28)] 6 . The 
Berry curvature \ is wirtten as d qa A\ — d\A qa . Under 
the periodic gauge, A\ is periodic in q a , and integration 
of d qa A\ over q a vanishes. Hence 



AP n 



BZ 



d( i , n 1 

(27r) d * 



A=0 



(2.28) 



In view of Eq. (2.28), both the adiabatic transport and 
the electric polarization can be regarded as the manifes- 
tation of Zak's phase, given in Eq. (1.29). 

However, a price must be paid to use the two-point for- 
mula, namely, the polarization in Eq. (2.28) is determined 
up to an uncertainty quantum. Since the integral (2.28) 
does not track the history of A, there is no information on 
how many cycles A has gone through. According to our 
discussion on quantized particle transport in Sec. II. B, for 
each cycle an integer number of electrons are transported 
across the sample, hence the polarization is changed by 
multiple of the quantum 



ea 



(2.29) 



where a is the Bravais lattice vector and Vo is the volume 
of the unit cell. 

Because of this uncertainty quantum, the polarization 
may be regarded as a multi-valued quantity with each 
value differed by the quantum. With this in mind, let us 
consider Zak's phase in a one-dimensional system with 
inversion symmetry. Now we know that Zak's phase is 
just 27r/e times the polarization density P. Under spa- 
tial inversion, P transforms to —P. On the other hand, 
inversion symmetry requires that P and — P describes 
the same state, which is only possible if P and — P differ 
by multiple of the polarization quantum ea. Therefore 



6 A more general phase choice is given by the path-independent 
gauge |n n (q,A)> = e^ e ^+ G ^\u n (q + G, X)) , where 6{q) is an 
arbitrary phase (Ortiz and Martin, 1994) 



drr\W n {r)f 



(2.30) 



where W n (r) is the Wannier function of the nth band, 
W n (r -R) = VNVo [ jP^e^-^Unkir) . (2.31) 

JBZ ( 2 V 

In this definition, one effectively maps a band insulator 
into a periodic array of localized distributions with truly 
quantized charges. This resembles an ideal ionic crystal 
where the polarization can be understood in the classical 
picture of localized charges. The quantum uncertainty 
found in Eq. (2.29) is reflected by the fact that the Wan- 
nier center position is defined only up to a lattice vector. 

Before concluding, we point out that the polar- 
ization defined above is clearly a bulk quantity as 
it is given by the Berry phase of the ground state 
wave function. A many-body formulation was devel- 
oped by Ortiz and Martin (1994) based on the work of 
Niu and Thouless (1984). 

Recent development in this field falls into two cate- 
gories. On the computational side, calculating polariza- 
tion in finite electric fields has been addressed, which has 
a deep influence on density functional theory in extended 
systems (Nunes and Gonze, 2001; Nunes and Vandcrbilt, 
1994; Souza et ai, 2002). On the theory side, Resta 
(1998) proposed a quantum-mechanical position operator 
for extended systems. It was shown that the expectation 
value of such an operator can be used to characterize 
the phase transition between the metallic and insulating 
state (Resta and Sorella, 1999; Souza et ai, 2000), and 
is closely related to the phenomenon of electron localiza- 
tion (Kohn, 1964). 

1. The Rice-Mele model 

So far our discussion of the adiabatic transport and 
electric polarization has been rather abstract. We now 
consider a concrete example: a one-dimensional dimcr- 
ized lattice model described by the following Hamiltonian 

H = + (-^'^Vi + h ' c -) + A(-l)'c}c,- +1 , 

3 

(2.32) 

where t is the uniform hopping amplitude, 5 is the 
dimerization order, and A is a staggered sublattice po- 
tential. It is the prototype Hamiltonian for a class 
of one-dimensional ferroelectrics. At half-filling, the 
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FIG. 2 (color online). Polarization as a function of A and 
S in the Rice-Mele model. The units is ea with a being the 
lattice constant. 



system is a metal for A = 5 = 0, and an insula- 
tor otherwise. Rice and Mele (1982) considered this 
model in the study of solitons in polyenes. It was 
later used to study ferroelectricity (Onoda et ai, 2004b; 
Vandcrbilt and King-Smith, 1993). If A = it reduces 
to the celebrated Su-Shrieffer-Heeger model (Su et ai, 
1979). 

Assuming periodic boundary condition, the Bloch rep- 
resentation of the above Hamiltonian is given by H(q) — 
h(q) ■ er, where 

h= (tcosy,-5sin^,A) . (2.33) 

This is the two- level model we discussed in Sec. I.C.3. Its 
energy spectrum consists of two bands with eigenenergies 
e± = ±(A 2 + S 2 sin 2 f + i 2 cos 2 f fl 2 . The degeneracy 
point occurs at 



• , (5 = 0, q = ir/a . 



(2.34) 



The polarization is calculated using the two-point for- 
mula (2.28) with the Berry connection given by 



Aq Oq 



dqOAe = sin 2 -d q (, 



(2.35) 



where 9 and <fr are the spherical angles of h. 

Let us consider the case of A = 0. In the parameter 
space of h, it lies in the xy-plane with 9 = tt/2. As q 
varies from to 2ir/a, tj> changes from to ir if 6 < 
and to — ir if S > 0. Therefore the polarization differ- 
ence between P(6) and P(—S) is ea/2. This is consistent 
with the observation that the state with P(—S) can be 
obtained by shifting the state with P(8) by half of the 
unit cell length a. 

Figure 2 shows the calculated polarization for arbi- 
trary A and S. As we can see, if the system adiabatically 
evolves along a loop enclosing the degeneracy point (0, 0) 



in the (A, 5) space, then the polarization will be changed 
by ea, which means that if we allow (A, 5) to change in 
time along this loop, for example, A(i) = A sin(i) and 
5(t) = Sq cos(t), a quantized charge of e is pumped out of 
the system after one cycle. On the other hand, if the loop 
does not contain the degeneracy point, then the pumped 
charge is zero. 



III. ELECTRON DYNAMICS IN THE PRESENCE OF 
ELECTRIC FIELDS 

The dynamics of Bloch electrons under the perturba- 
tion of an electric field is one of the oldest problems in 
solid state physics. It is usually understood that while 
the electric field can drive electron motion in the mo- 
mentum space, it does not appear in the electron ve- 
locity; the latter is simply given by (for example, see 
Ashcroft and Mcrmin, 1976) 



v n (q) 



de n (q) 
fidq 



(3.1) 



Through recent progress on the semiclassical dynamics 
of Bloch electrons, it has been made increasingly clear 
that this description is incomplete. In the presence 
of an electric field, an electron can acquire an anoma- 
lous velocity proportional to the Berry curvature of the 
band (Chang and Niu, 1995, 1996; Sundaram and Niu, 
1999). This anomalous velocity is responsible for a num- 
ber of transport phenomena, in particular various Hall 
effects, which we study in this section. 



A. Anomalous velocity 

Let us consider a crystal under the perturbation of a 
weak electric field E, which enters into the Hamiltonian 
through the coupling to the electrostatic potential 4>(r). 
However, a uniform E means that 4>{r) varies linearly 
in space and breaks the translational symmetry of the 
crystal such that Bloch's theorem cannot be applied. To 
go around this difficulty, one can let the electric field enter 
through a uniform vector potential A(t) that changes in 
time. Using the Peierls substitution, the Hamiltonian is 
written as 



H(t) 



[P + eA{t)f 
2m 



V(r) 



(3.2) 



This is the time-dependent problem we have studied in 
last section. Transforming to the q-space representation, 
we have 

H(q,t) = H(q+^A(t)) . (3.3) 
Introduce the gauge-invariant crystal momentum 

k = q+^A(t). (3.4) 
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The parameter-dependent Hamiltonian can be simply 
written as H(k(q,t)). Hence the eigenstates of the time- 
dependent Hamiltonian can be labeled by a single pa- 
rameter fe. Moreover, because A(t) preserves the transla- 
tional symmetry, q is still a good quantum number and is 
a constant of motion q — 0. It then follows from Eq. (3.4) 
that k satisfies the following equation of motion 

k = ~E. (3.5) 

Using the relation d/dq a = d/dk a and d/dt = 
— (e/Ti)E a d/dk a , the general formula (2.5) for the ve- 
locity in a given state k becomes 

v n (k) = %P - € -E x n B (fc) , (3.6) 

where fl n (k) is the Berry curvature of the nth band: 

n n {k) = i{v (k)\ x |V fc u„(fe)> . (3.7) 

We can see that, in addition to the usual band dis- 
persion contribution, an extra term previously known 
as an anomalous velocity also contributes to v n (k). 
This velocity is always transverse to the electric 
field, which will give rise to a Hall current. His- 
torically, the anomalous velocity has been obtained 
by Adams and Blount (1959); Karplus and Luttingcr 
(1954); Kohn and Luttinger (1957); however, its relation 
to the Berry phase was not realized. In Sec. V we shall 
rederive this term using a wave packet approach. 

B. Berry curvature: Symmetry considerations 

The velocity formula (3.6) reveals that, in addition 
to the band energy, the Berry curvature of the Bloch 
bands is also required for a complete description of the 
electron dynamics. However, the conventional formula, 
Eq. (3.1), has had great success in describing various elec- 
tronic properties in the past. It is thus important to know 
under what conditions the Berry curvature term cannot 
be neglected. 

The general form of the Berry curvature i~2 n (fc) can 
be obtained via symmetry analysis. The velocity for- 
mula (3.6) should be invariant under time reversal and 
spatial inversion operations if the unperturbed system 
has these symmetries. Under time reversal, v n and k 
changes sign while E is fixed. Under spatial inversion, 
v n , k, and E changes sign. If the system has time re- 
versal symmetry, the symmetry condition on Eq. (3.6) 
requires that 

n„(-fc) = -n„(fc) . (3.8) 

If the system has spatial inversion symmetry, then 

n„(-fe) = n„(fe) . (3.9) 

Therefore, for crystals with simultaneous time-reversal 
and spatial inversion symmetry the Berry curvature van- 
ish identically throughout the Brillouin zone. In this case 
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FIG. 3 (color online). Fermi surface in (010) plane (solid 
lines) and Berry curvature — Sl z (k) in atomic units (color 
map) of fee Fe. From Yao et ai, 2004. 

Eq. (3.6) reduces to the simple expression (3.1). However, 
in systems with either broken time-reversal or inversion 
symmetries, their proper description requires the use of 
the full velocity formula (3.6). 

There are many important physical systems where 
both symmetries are not simultaneously present. For 
example, in the presence of ferro- or antiferro-magnetic 
ordering the crystal breaks the time-reversal symmetry. 
Figure 3 shows the Berry curvature on the Fermi sur- 
face of fee Fe. As we can see, the Berry curvature is 
negligible in most areas in the momentum space and dis- 
plays very sharp and pronounced peaks in regions where 
the Fermi lines (intersection of the Fermi surface with 
(010) plane) have avoided crossings due to spin-orbit cou- 
pling. The monopolc structure has been identified in 
other materials as well (Fang et ai, 2003). Another ex- 
ample is provided by single-layered graphene sheet with 
staggered sublattice potential, which breaks the inver- 
sion symmetry (Zhou et ai, 2007). Figure 4 shows the 
energy band and Berry curvature of this system. The 
Berry curvature at valley K\ and if 2 have opposite signs 
due to time-reversal symmetry. We note that as the gap 
approaches zero, the Berry phase acquired by an elec- 
tron during one circle around the valley becomes exactly 
±7r. This Berry phase of 7r has been observed in intrin- 
sic graphene sheet (Novosclov et ai, 2005; Zhang et ai, 
2005). 

C. The quantum Hall efFect 

The quantum Hall effect was discovered by 
Klitzing et al. (1980). They found that in a strong 
magnetic field the Hall conductivity of a two-dimensional 
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FIG. 4 (color online). Energy bands (top panel) and 
Berry curvature of the conduction band (bottom panel) of 
a graphene sheet with broken inversion symmetry. The first 
Brillouin zone is outlined by the dashed lines, and two inequiv- 
alent valleys are labeled as Ki and K2. Details are presented 
in Xiao et al, 2007b. 



electron gas is exactly quantized in the units of e 2 /h. 
The exact quantization was subsequently explained 
by Laughlin (1981) based on gauge invariance, and 
was later related to a topological invariance of the 
energy bands (Avron et al., 1983; Niu et al., 1985; 
Thouless et al., 1982). Since then it has blossomed into 
an important research field in condensed matter physics. 
In this subsection we shall focus only on the quantization 
aspect of the quantum Hall effect using the formulation 
developed so far. 

Let us consider a two-dimensional band insulator. It 
follows from Eq. (3.6) that the Hall conductivity of the 
system is given by 



>xy 



BZ 



d 2 k 
(2^ 



(3.10) 



where the integration is over the entire Brillouin. Once 
again we encounter the situation where the Berry curva- 
ture is integrated over a closed manifold. Here a xy is just 
the Chern number in the units of e 2 /h, i.e., 

e 2 

<?x V = n— . (3-11) 

Therefore the Hall conductivity is quantized for a two- 
dimensional band insulator of non-interacting electrons. 

Historically, the quantization of the Hall conductivity 
in a crystal was first shown by Thouless et al. (1982) for 
magnetic Bloch bands (see also Sec. VIII. It was shown 
that, due to the magnetic translational symmetry, the 
phase of the wave function in the magnetic Brillouin zone 
carries a vortex and leads to a nonzero quantized Hall 
conductivity (Kohmoto, 1985). However, from the above 
derivation it is clear that a magnetic field is not necessary 



for the quantum Hall effect to occur as the condition is 
a nonzero Chern number of the band. Haldane (1988) 
constructed a tight-binding model on a honeycomb lat- 
tice which displays the quantum Hall effect with zero 
net flux per unit cell. Another model utilizing the spin- 
orbit interaction in a semiconductor quantum well was 
recently proposed (Liu et al., 2008; Qi et al., 2006). The 
possibility of realizing the quantum Hall effect without a 
magnetic field is very attractive in device design. 

Niu et al. (1985) further showed that the quantized 
Hall conductivity in two-dimensions is robust against 
many-body interactions and disorder. Their derivation 
involves the same technique discussed in Sec. II. B. 2. A 
two-dimensional many-body system is placed on a torus 
by assuming periodic boundary conditions in both di- 
rections. One can then thread the torus with magnetic 
flux through its holes (Fig. 5) and make the Hamilto- 
nian H{4>It4>2) depend on the flux tfii and <f>2- The Hall 
conductivity is calculated using the Kubo formula 
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(3.12) 

where is the many-body wave function with |$q) the 
ground state. In the presence of the flux, the veloc- 
ity operator is given by Uj = dH (k 1; K^j/difiKi) with 
Ki = (e/K)(f>i/Li and Li the dimensions of the sys- 
tem. We recognize that Eq. (3.12) is the summation for- 
mula (1.13) for the Berry curvature Q KlK2 of the state 
|<f> ). The existence of a bulk energy gap guarantees that 
the Hall conductivity remains unchanged after thermo- 
dynamic averaging, which is given by 



a H = 



2-n/Lt 



dux 



2tt/L 2 



dn2 fi„ 



(3.13) 



Note that the Hamiltonian H(k\,K2) is periodic in m 
with period 2n/Li because the system returns to its orig- 
inal state after the flux is changed by a flux quantum h/e 
(and Ki changed by 2w/Li). Therefore the Hall conduc- 
tivity is quantized even in the presence of many-body 
interaction and disorder. Due to the high precision of 
the measurement and the robustness of the quantization, 
the quantum Hall resistance is now used as the primary 
standard of resistance. 

The geometric and topological ideas developed in the 
study of the quantum Hall effect has a far-reaching im- 
pact on modern condensed matter physics. The robust- 
ness of the Hall conductivity suggests that it can be used 
as a topological invariance to classify many-body phases 
of electronic states with a bulk energy gap (Avron et al. , 
1983): states with different topological orders (Hall con- 
ductivities in the quantum Hall effect) cannot be adia- 
batically transformed into each other; if that happens, a 
phase transition must occur. It has important applica- 
tions in strongly correlated electron systems, such as the 
fractional quantum Hall effect (Wen and Niu, 1990), and 
most recently, the topological quantum computing (for a 
review, see Nayak et al., 2008). 
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FIG. 5 Magnetic flux going through the holes of the torus. 



D. The anomalous Hall Effect 

Next we discuss the anomalous Hall effect, which 
refers to the appearance of a large spontaneous Hall cur- 
rent in a ferromagnet in response to an electric field 
alone (Chien and Westgate, 1980). Despite its century- 
long history and importance in sample characterization, 
the microscopic mechanism of the anomalous Hall effect 
has been a controversial subject (for a brief review on 
the history of the anomalous Hall effect, see Sinova et ai, 
2004). In the past, three mechanisms have been identi- 
fied: the intrinsic contribution (Karplus and Luttinger, 
1954; Luttinger, 1958), and the extrinsic contributions 
from the skew (Smit, 1958) and side-jump scatter- 
ing (Berger, 1970). It was later realized that the 
scattering-independent intrinsic contribution comes from 
the Berry-phase supported anomalous velocity. It will be 
our primary interests here. 

The intrinsic contribution to the anomalous Hall ef- 
fect can be regarded as an "unquantized" version of the 
quantum Hall effect. The Hall conductivity is given by 



dk 

h J {2nj 



(3.14) 



where /(£&) is the Fermi-Dirac distribution function. 
However, unlike the quantum Hall effect, the anomalous 
Hall effect does not require a nonzero Chern number of 
the band; for a band with zero Chern number, the local 
Berry curvature can be nonzero and give rise to a nonzero 
anomalous Hall conductivity. 

To appreciate the intrinsic contribution, let us con- 
sider a generic Hamiltonian with spin-orbit (SO) split 
bands (Onoda et ai, 2006b) 

Ti 2 k 2 

H= h A(fe x er) • e z - Aa z . (3.15) 

2m 

where 2A is the SO split gap in the energy spectrum 
e± = h 2 k 2 /2m ± \/A 2 fc 2 + A 2 , and A gives a linear dis- 
persion in the absence of A. This model also describes 
spin-polarized two-dimensional electron gas with Rashba 
SO coupling, with A being the SO coupling strength and 
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FIG. 6 (a) Energy dispersion of spin-split bands, (b) The 
Hall conductivity — a xy in the units of e 2 /h as a function of 
Fermi energy. 



A the exchange field (Culcer et al., 2003). 
curvature is given by, using Eq. (1.19), 



n± = T 



A 2 A 



2(A 2 fc 2 + A 2 ) 3 / 2 



The Berry 



(3.16) 



The Berry curvatures of the two energy bands have oppo- 
site sign, and is highly concentrated around the gap (In 
fact, the Berry curvature has the same form of the Berry 
curvature in one valley of the graphene, shown in Fig. 4). 
One can verify that the integration of the Berry curva- 
ture of a full band, 27r qdqfl±, is ±7r for the upper 
and lower bands, respectively. 

Figure 6 shows the band dispersion, and the intrin- 
sic Hall conductivity, Eq. (3.14), as the Fermi energy 
sweeps across the SO split gap. As we can see, when 
the Fermi energy sf is in the gap region, the Hall con- 
ductivity reaches its maximum value (about —e 2 /2h). If 
£f < —A, the states with energies just below —A, which 
contribute most to the Hall conductivity, are empty. If 
£p > A, contributions from upper and lower bands cancel 
each other, and the Hall conductivity decreases quickly 
as £f moves away from the bad gap. It is only when 
— A < Ep < A, the Hall conductivity is resonantly en- 
hanced (Onoda et ai, 2006b). 
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1. Intrinsic vs. extrinsic contributions 

The above discussion does not take into account the 
fact that, unlike insulators, in metallic systems electron 
scattering can be important in transport phenomena. 
Two contributions to the anomalous Hall effect arises 
due to scattering: (i) the skew scattering that refers to 
the asymmetric scattering amplitude with respect to the 
scattering angle between the incoming and outgoing elec- 
tron waves (Smit, 1958), and (ii) the side jump which is 
a sudden shift of the electron coordinates during scat- 
tering (Berger, 1970). In a more careful analysis, a sys- 
tematic study of the anomalous Hall effect based on the 
semiclassical Boltzmann transport theory has been car- 
ried out (Sinitsyn, 2008). The basic idea is to solve the 
following Boltzmann equation: 



dgk 
dt 



eE 



de df 
fidk de 



E 



9k ~ 9k' 



-%eE- Hi- 
de 



hh 1 



(3.17) 



where g is the non-equilibrium distribution function, ujkk' 
represents the asymmetric skew scattering, and Sr^k' de- 
scribes the side-jump of the scattered electrons. The Hall 
conductivity is the sum of different contributions 



oh = cr H 



' H 



'H > 



(3.18) 



where is the intrinsic contribution given by Eq. (3.14), 
crff is the skew scattering contribution, which is propor- 
tional to the relaxation time t, and a H is the side jump 
contribution, which is independent of r. Note that in 
addition to Berger 's original proposal, cr^ also includes 
two other contributions: the intrinsic skew-scattering and 
anomalous distribution function (Sinitsyn, 2008). 

An important question is to identify the dominant con- 
tribution to the AHE among these mechanisms. If the 
sample is clean and the temperature is low, the relaxation 
time r can be extremely large, and the skew scattering is 
expected to dominate. On the other hand, in dirty sam- 
ples and at high temperatures, becomes small com- 
pared to both a H and a H . Because the Berry-phase con- 
tribution a h is independent of scattering, it can be read- 
ily evaluated using first-principles methods or effective 
Hamiltonians. Excellent agreement with experiments has 
been demonstrated in ferromagnetic transition metals 
and semiconductors (Fang et al, 2003; Jungwirth et al, 
2002; Xiao et al, 2006b; Yao et al, 2004, 2007), which 
leaves little room for the side jump contribution. 

In addition, a number of experimental results also 
gave favorable evidence for the dominance of the in- 
trinsic contribution (Chun et al, 2007; Lee et al, 2004b; 
Mathieu et al, 2004; Sales et al, 2006; Zeng et al, 
2006). In particular, Tian et al. recently measured the 
anomalous Hall conductivity in Fe thin films. By varying 
the film thickness and the temperature, they are able to 
control various scattering process such as the impurity 
scattering and the phonon scattering. Figure 7 shows 
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FIG. 7 (color online) a a h vs. cr xx (T) 2 in Fe thin films with 
different film thickness over the temperature range 5-300 K. 
From Tian et al. 



their measured a a h as a function of a xx (T) 2 . We can see 
that although a a h in different thin films and at different 
temperatures shows a large variation at finite a xx , they 
converge to a single point as a xx approaches zero, where 
the impurity-scattering induced contribution should be 
washed out by the phonon scattering and only the intrin- 
sic contribution survives. It turns out that this converged 
value is very close to the bulk a a h of Fe, which confirms 
the dominance of the intrinsic contribution in Fe. 

In addition to the semiclassical approach (Sinitsyn, 
2008; Sinitsyn et al, 2005), there are a number 
of works based on a full quantum mechanical ap- 
proach (Dugaev et al, 2005; Inoue et al, 2006 
Kato et al, 2007; Nozieres and Lewiner, 1973 
Onoda and Nagaosa, 2002; Onoda et al, 2006b, 2008 
Sinitsyn et al, 2007). In both approaches, the Berry- 
phase supported intrinsic contribution to the anomalous 
Hall effect has been firmly established. 



2. Anomalous Hall conductivity as a Fermi surface property 

An interesting aspect of the intrinsic contribution to 
the anomalous Hall effect is that the Hall conductivity, 
Eq. (3.14), is given as an integration over all occupied 
states below the Fermi energy. It seems to be against the 
spirit of the Landau Fermi liquid theory, which states 
that the transport property of an electron system is de- 
termined by quasiparticles at the Fermi energy. This 
issue was first raised by Haldane (2004), and he showed 
that the Hall conductivity can be written, in the units of 
e 2 /2irh, as the Berry phase of quasiparticles on the Fermi 
surface, up to a multiple of 2tt. Therefore the intrinsic 
Hall conductivity is also a Fermi surface property. This 
observation suggests that the Berry phase on the Fermi 
surface should be added as a topological ingredient to the 
Landau Fermi liquid theory. 

For simplicity, let us consider a two-dimensional sys- 
tem. We assume there is only one partially filled band. 
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Write the Berry curvature in terms of the Berry vector 
potential and integrate Eq. (3.14) by part; one finds 

Note that the Fermi distribution function / is a step func- 
tion at the Fermi energy. If we assume the Fermi surface 
is a closed loop in the Brillouin zone, then 

*ly = ^ j> dk ■ A k . (3.20) 

The integral is nothing but the Berry phase along the 
Fermi circle in the Brillouin zone. The three-dimensional 
case is more complicated; Haldane (2004) showed that 
the same conclusion can be reached. 

Wang et al. (2007) has implemented Haldane 's idea in 
first-principles calculations of the anomalous Hall con- 
ductivity. From a computational point of view, the ad- 
vantage lies in that the integral over the Fermi sea is con- 
verted to a more efficient integral on the Fermi surface. 
On the theory side, Shindou and Balents (2006, 2008) 
derived an effective Boltzmann equation for quasiparti- 
cles on the Fermi surface using the Keldysh formalism, 
where the Berry phase of the Fermi surface is defined in 
terms of the quasiparticle Green functions, which nicely 
fits into the Landau Fermi liquid theory. 

E. The valley Hall efFect 

A necessary condition for the charge Hall effect to man- 
ifest is the broken time-reversal symmetry of the system. 
In this subsection we discuss another type of Hall effect 
which relies on inversion symmetry breaking, and sur- 
vives in time-reversal invariant systems. 

We shall use graphene as our prototype system. 
The band structure of intrinsic graphene has two de- 
generate and inequivalcnt Dirac points at the cor- 
ners of the Brillouin zone, where the conduction and 
valance bands touch each other, forming a valley struc- 
ture. Because of their large separation in momen- 
tum space, the intervalley scattering is strongly sup- 
pressed (Gorbachev et al, 2007; Morozov et al, 2006; 
Morpurgo and Guinea, 2006), which makes the valley 
index a good quantum number. Interesting valley- 
dependent phenomena, which concerns about the de- 
tection and generation of valley polarization, are be- 
ing actively explored (Akhmcrov and Bccnakkcr, 2007; 
Rycerz et al, 2007; Xiao et al, 2007b; Yao et al, 2008). 

The system we are interested in is graphene with bro- 
ken inversion symmetry. Zhou et al. (2007) have re- 
cently reported the observation of a band gap open- 
ing in epitaxial graphene, attributed to the inversion 
symmetry breaking by the substrate potential. In 
addition, in biased graphene bilayer, inversion sym- 
metry can be explicitly broken by the applied in- 
terlayer voltage (McCann and Fal'ko, 2006; Min et al, 



2007; Ohta et al., 2006). It is this broken inversion sym- 
metry that allows a valley Hall effect. Introducing the 
valley index t z — ±1 which labels the two valleys, we 
can write the valley Hall effect as 

f = <t v h zxE, (3.21) 

where a v H is the valley Hall conductivity, and the valley 
current j v = (r z v) is defined as the average of the val- 
ley index t z times the velocity operator v. Under time 
reversal, both the valley current and electric field are in- 
variant (t z changes sign because the two valleys switch 
when the crystal momentum changes sign). Under spa- 
tial inversion, the valley current is still invariant but the 
electric field changes sign. Therefore, the valley Hall con- 
ductivity can be nonzero when the inversion symmetry 
is broken, even if time reversal symmetry remains. 

In the tight-binding approximation, the Hamiltonian of 
a single graphene sheet can be modeled with a nearest- 
neighbor hopping energy t and a site energy difference A 
between sublattices. For relatively low doping, we can re- 
sort to the low-energy description near the Dirac points. 
The Hamiltonian is given by (Semenoff, 1984) 

\/3 A 
H = —at(q x T z a x + q y a y ) + —a z , (3.22) 

where <x is the Pauli matrix accounting for the sublattice 
index, and q is measured from the valley center If 1.2 = 
(±47r/3a)a; with a being the lattice constant. The Berry 
curvature of the conduction band is given by 

Note that the Berry curvatures in two valleys have oppo- 
site sign, as required by time-reversal symmetry. The en- 
ergy spectrum and the Berry curvature are already shown 
in Fig. 4. Once the structure of the Berry curvature is re- 
vealed, the valley Hall effect becomes transparent: upon 
the application of an electric field, electrons in different 
valleys will flow to opposite transverse edges, giving rise 
to a net valley Hall current in the bulk. 

We remark that as A goes to zero, the Berry curvature 
vanishes everywhere except at the Dirac points where it 
diverges. Meanwhile, The Berry phase around the Dirac 
points becomes exactly ±7r (also see Sec. VII. C). 

As we can see, the valley transport in systems 
with broken inversion symmetry is a very general phe- 
nomenon. It provides a new and standard pathway 
to potential applications of valley tronics, or valley- 
based electronic applications, in a broad class of semi- 
conductors (Gunawan et al., 2006; Xiao et al, 2007b; 
Yao et al, 2008). 

IV. WAVE PACKET: CONSTRUCTION AND 
PROPERTIES 

So far, our discussion has focused on crystals under 
time-dependent perturbations, and we have shown that 
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the Berry phase will manifest itself as an anomalous term 
in the electron velocity. However, in general situations 
the electron dynamics can be also driven by perturba- 
tions that vary in space. In this case, the most use- 
ful description is provided by the semiclassical theory 
of Bloch electron dynamics, which describes the motion 
of a narrow wave packet obtained by superposing the 
Bloch states of a band (see, for example, Chap. 12 of 
Ashcroft and Mermin, 1976). The current and next sec- 
tions are devoted to the study of the wave packet dynam- 
ics, where the Berry curvature naturally appears in the 
equations of motion. 

In this section we discuss the construction and the 
general properties of the wave packet. Two quantities 
emerges in the wave packet approach, i.e., the orbital 
magnetic moment of the wave packet and the dipole mo- 
ment of a physical observable. For their applications, 
we consider the problems of orbital magnetization and 
anomalous thermoelectric transport in ferromagnets. 

A. Construction of the wave packet and its orbital moment 

We assume the perturbations are sufficiently weak such 
that transitions between different bands can be neglected, 
i.e., the electron dynamics is confined within a single 
band. Under this assumption, we construct a wave packet 
using the Bloch functions \ip n (q)) from the nth band: 

\W ) = J dqw(q,t)\Mq)) ■ (4-1) 

There are two requirements on the envelope function 
w(q,t). Firstly, w(q,t) must have a sharp distribution 
in the Brillouin zone such that it makes sense to speak 
of the wave vector q c of the wave packet, given by 

Qc = J dqq\w(q,t)\ 2 . (4.2) 

To first order, \w(q, t)| 2 can be approximated by 5(q — q c ) 
and one has 

Jdqf(q)\w( qi t)\ 2 ^f(q c ), (4.3) 

where f(q) is an arbitrary function of q. Equation (4.3) 
is very useful in evaluating various quantities related to 
the wave packet. Secondly, the wave packet has to be 
narrowly localized around its center of mass, denoted by 
r c , in the real space, i.e., 

r c = (W Q \r\W ) . (4.4) 

Using Eq. (4.3) we obtain 

d 

r c = -g— axgw(q c ,t) + A£(q c ) , (4.5) 

where A!q — i(u n (q)\ V q |n„(q)) is the Berry connection 
of the nth band defined using \u n (q)) — e~ r \ip n (q)). 



In principle, more dynamical variables, such as the 
width of the wave packet in both the real space and 
momentum space, can be added to allow a more elabo- 
rate description of the time evolution of the wave packet. 
However, in short period the dynamics is dominated by 
the motion of the wave packet center, and Eqs. (4.2) and 
(4.5) are sufficient requirements. 

When more than one band come close to each other, or 
even become degenerate, the single-band approximation 
breaks down and the wave packet must be constructed 
using Bloch functions from multiple bands. Culccr et al. 
(2005); Shindou and Imura (2005) developed the multi- 
band formalism for electron dynamics, which will be pre- 
sented in Sec. IX. For now, we will focus on the single- 
band formulation and drop the band index n for simple 
notation. 

The wave packet, unlike a classical point particle, has 
a finite spread in real space. In fact, since it is con- 
structed using an incomplete basis of the Bloch func- 
tions, the size of the wave packet has a nonzero lower 
bound (Marzari and Vandcrbilt, 1997). Therefore, a 
wave packet may possess a self-rotation around its center 
of mass, which will in turn give rise to an orbital mag- 
netic moment. By calculating the angular momentum of 
a wave packet directly, one finds (Chang and Niu, 1996) 

rn(q) = -*(W \(r-r c ) x j\W ) 

% (4-6) 
= - l ^( y i u \ x [#(<?) " e(*l)]|V,ti} , 

where H(q) = e~ lq r H e lq r is the q-dependcnt Hamilto- 
nian. This shows that the wave packet of a Bloch electron 
generally rotates around its mass center and carries an 
orbital magnetic moment in addition to its spin moment. 

We emphasize that the orbital moment is an intrin- 
sic property of the band. Its final expression, Eq. (4.6), 
does not depend on the actual shape and size of the wave 
packet, and only depends on the Bloch functions. Un- 
der symmetry operations, the orbital moment transforms 
exactly like the Berry curvature. Therefore unless the 
system has both time-reversal and inversion symmetry, 
m(q) is in general nonzero. Information of the orbital 
moment can be obtained by measuring magnetic circular 
dichroism spectrum of a crystal (Souza and Vandcrbilt, 
2008; Yao et al, 2008). 

This orbital moment behaves exactly like the electron 
spin. For example, upon the application of a magnetic 
field, the orbital moment will couple to the field through 
a Zeeman-like term — m(q) ■ B. If one construct a wave 
packet using only the positive energy states (the electron 
branch) of the Dirac Hamiltonian, its orbital moment in 
the non-relativistic limit is exactly the Bohr magneton 
(Sec. IX). For Bloch electrons, the orbital moment can 
be related to the electron factor (Yafet, 1963). Let 
us consider a specific example. For the graphene model 
with broken-inversion symmetry, discussed in Sec. III.E, 
the orbital moment of the conduction band is given 
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by (Xiao et ai, 2007b) 

. . 3ea 2 At 2 fA . 

So orbital moments in different valleys have opposite 
signs, as required by time-reversal symmetry. Interest- 
ingly, the orbital moment at exactly the band bottom 
takes the following form 

eh 

m(T z ) = T z fj,* B , fi* B = ^—^ , (4.8) 

where m* is the effective mass at the band bottom. The 
close analogy with the Bohr magneton for the electron 
spin is transparent. In realistic situations, the moment 
can be 30 times larger than the spin moment, and can be 
used as an effective way to detect and generate the valley 
polarization (Xiao et ai, 2007b; Yao et ai, 2008). 

B. Orbital magnetization 

A closely related quantity to the orbital magnetic mo- 
ment is the orbital magnetization in a crystal. Al- 
though this phenomenon has been known for a long time, 
our understanding of orbital magnetization in crystals 
has remained in a primitive stage. In fact, there was 
no proper way to calculate this quantity until recently 
when the Berry phase theory of orbital magnetization 
is developed (Shi et ai, 2007; Thonhauser et ah, 2005; 
Xiao et al., 2005). Here we provide a rather pictorial 
derivation of the orbital magnetization based on the wave 
packet approach. This derivation gives an intuitive pic- 
ture of different contributions to the total orbital mag- 
netization. 

The main difficulty of calculating the orbital magneti- 
zation is exactly the same one we faced when calculating 
the electric polarization: the magnetic dipole er x p is 
not defined in a periodic system. For a wave packet this 
is not a problem because it is localized in space. As we 
showed in last subsection, each wave packet carries an or- 
bital moment. Thus, it is tempting to conclude that the 
orbital magnetization is simply the thermodynamic av- 
erage of the orbital moment. As it turns out, this is only 
part of the contribution. There is another contribution 
due to the center-of-mass motion of the wave packet. 

For simplicity, let us consider a finite system of elec- 
trons in a two-dimensional lattice with a confining po- 
tential V(r). We further assume that the potential V(r) 
varies slowly at atomic length scale such that the wave 
packet description of the electron is still valid on the 
boundary. In the bulk where V(r) vanishes, the electron 
energy is just the bulk band-energy; near the boundary, 
it will be tilted up due to the increase of V(r). Thus to 
a good approximation, we can write the electron energy 
as 

e(r,q) = s(q) + V(r). (4.9) 






x 

FIG. 8 Electron energy e in a slowly varying confining po- 
tential V(r). In addition to the self-rotation, wave packets 
near the boundary will also move along the boundary due 
to the potential V. Level spacings between different bulk q- 
states are exaggerated; they are continuous in the semiclassi- 
cal limit. The insert shows directions of the Berry curvature, 
the effective force, and the current carried by a wave packet 
on the left boundary. 

The energy spectrum in real space is sketched in Fig. 8. 

Before proceeding further, we need generalize the ve- 
locity formula (3.6), which is derived in the presence of 
an electric field. In our derivation the electric field en- 
ters through a time-dependent vector potential A(t) so 
that we can avoid the technical difficulty of calculating 
the matrix element of the position operator. However, 
the electric field may be also given by the gradient of the 
electrostatic potential. In both cases, the velocity for- 
mula should stay the same because it should be gauge 
invariant. Therefore, in general a scalar potential V(r) 
will induce a transverse velocity of the following form 

-W(r) x n(q) . (4.10) 

This generalization will be justified in Sec. VI 

Now consider a wave packet in the boundary region. It 
will feel a force Vl'(r) due to the presence of the confin- 
ing potential. Consequently, according to Eq. (4.10) the 
electron acquires a transverse velocity, whose direction is 
parallel with the boundary (Fig. 8). This transverse ve- 
locity will lead to a boundary current (of the dimension 
"current density x width" in 2D) given by 

where x is in the direction perpendicular to the boundary, 
and the integration range is taken from deep into the bulk 
to outside the sample. Recall that for a current / flowing 
in a closed circuit enclosing a sufficiently small area A, 
the circuit carries a magnetic moment given by I ■ A. 
Therefore the magnetization (magnetic moment per unit 
area) has the magnitude of the current /. Integrating 
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Eq. (4.11) by part, we obtain 

M f = - e Jdef(e)a xy (e), (4.12) 

where a xy (e) is the zero-temperature Hall conductivity 
for a system with Fermi energy e: 

a xy (e) = ~J ^©(e - e(q))n z (q) . (4.13) 

Since the boundary current corresponds to the global 
movement of the wave packet center, we call this contri- 
bution the "free current" contribution, whereas the or- 
bital moment are due to "localized" current. The total 
magnetization thus is 

M z = J j^y i f(q)m z (q)+ 1 - J ds f(e)a xy (e) . (4.14) 

The orbital magnetization has two different contribu- 
tions: one is from the self-rotation of the wave packet, 
and the other is due to the center-of-mass motion. 

The reader may still have one question in mind: the 
above derivation relies on the existence of a confining 
potential, which seems to contradict the fact that the or- 
bital magnetization is a bulk property. This is a wrong 
assertion as the final expression, Eq. (4.14), is given 
in terms of the bulk Bloch functions and does not de- 
pend on the boundary condition. Here, the boundary 
is merely a tool to expose the "free current" contribu- 
tion because in a uniform system, the magnetization 
current always vanishes in the bulk. Finally, in more 
rigorous approaches (Shi et ai, 2007; Xiao et ai, 2005) 
the boundary is not needed and the derivation is based 
on a pure bulk picture. It is similar to the quantum 
Hall effect, which can be understood in terms of ei- 
ther the bulk states (Thoulcss et ai, 1982) or the edge 
states (Halpcrin, 1982). 

C. Dipole moment 

The finite size of the wave packet not only allows an 
orbital magnetic moment, but also leads to the concept 
of the dipole moment associated with an operator. 

The dipole moment appears naturally when we con- 
sider the thermodynamic average of a physical quantity, 
with its operator denoted by O. In the wave packet ap- 
proach, it is given by 

0(r) = J ^0g(r c ,q c )(W\dS(r-r)\W) , (4.15) 

where g(r,q) is the distribution function, (W|««-|W) 
denotes the expectation in the wave packet state, and 
<5(r — r) plays the role as a sampling function, as shown in 
Fig. 9. An intuitive way to view Eq. (4.15) is to think the 
wave packets as small molecules, then Eq. (4.15) is the 
quantum mechanical version of the familiar coarse grain- 
ing process which averages over the length scale larger 



than the size of the wave packet. A multi-pole expansion 
can be carried out. But for most purposes, the dipole 
term is enough. Expand the (^-function to first order of 
r - r c : 

6(r - r) = S(r - r c ) - (r - r c ) ■ VS(r - r c ) . (4.16) 
Inserting it into Eq. (4.15) yields 

0(r)= J JjLg( r ,q)(W\d\W)\ Pe=r 

-V- J ^-g(r,q){W\d(f-r c )\W)\ rc=r . 

(4.17) 

The first term is what one would obtain if the wave packet 
is treated as a point particle. The second term is due to 
the finite size of the wave packet. We can see that the 
bracket in the second integral has the form of a dipole of 
the operator O, defined by 

P 6 = (W\6(f - r c )\W) , (4.18) 

The dipole moment of an observable is a general conse- 
quence of the wave packet approach and must be included 
in the semiclassical theory. Its contribution appears only 
when the system is inhomogeneous. 
In particular, we find: 

1. If O = e, then P e = 0. This is consistent with the 
fact that the charge center coincides with the mass 
center of the electron. 

2. If O = v, one finds the expression for the local 
current: 

jL = J (2^ 5(r ' q) * + V X / (2^- 9(r ' Q)m{q) ■ 

(4.19) 

We will explain the meaning of local later. Inter- 
estingly, this is the second time we encounter the 
quantity m(q), but in an entirely different con- 
text. The physical meaning of the second term 




FIG. 9 Sampling function and a wave packet at r c . The 
width L of the sampling function is sufficiently small so that 
it can be treated as a 5-function at the macroscopic level and 
is sufficiently big so that it contains a large number of wave 
packets of width I inside its range. Eq. (4.15) is indeed a mi- 
croscopic average over the distance L around the point r. See 
Section 6.6 in Jackson (1998) for an analogy in macroscopic 
electromagnetism. 
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wave packet 
W{r;r c ,kc) 



FIG. 10 The wave packet description of a charge carrier whose 
center is (r c , Qc). A wave packet generally possesses two kinds 
of motion: the center of mass motion and the self-rotation 
around its center. From Xiao et al, 2006b. 



becomes transparent if we make reference to the 
self-rotation of the wave packet. The self-rotation 
can be thought as localized circuit. Therefore if the 
distribution is not uniform, those localized circuit 
will contribute to the local current j L (See Fig. 10). 

3. If O is the spin operator s, then Eq. (4.18) gives 
the spin dipole 



P s = {u\s(i-^-A q ) 



(4.20) 



It shows that in general the spin center and the 
mass center do not coincide, which is usually due 
to the spin-orbit interaction. The time derivative 
of the spin dipole contributes to the total spin cur- 
rent (Culcer et al, 2004). 



D. Anomalous thermoelectric transport 

As an application of the above concepts, we con- 
sider the problem of anomalous thermoelectric trans- 
port in ferromagnets, which refers to the appearance 
of a Hall current driven by statistical forces, such 
as the gradient of temperature and chemical poten- 
tial (Chicn and Wcstgate, 1980). Similar to the anoma- 
lous Hall effect, there are also intrinsic and extrinsic con- 
tributions, and we will focus on the former. 

A question immediately arises when one tries to for- 
mulate this problem. Recall that in the presence of an 
electric field, the electron acquires an anomalous velocity 
proportional to the Berry curvature, which gives rise to 
a Hall current. In this case, the driving force is of me- 
chanical nature: it exists on the microscopic level and 
can be described by a perturbation to the Hamiltonian 
for the carriers. On the other hand, transport can be 
also driven by statistical force. However, the statisti- 
cal force manifests on the macroscopic level and makes 
sense only through the statistical distribution of the car- 
riers. Since there is no force acting directly on individual 
particles, the obvious cause for the Berry phase assisted 
transport is eliminated. This conclusion would introduce 
a number of basis contradictions to the standard trans- 
port theory. First, a chemical potential gradient would 
be distinct from the electric force, violating the basis for 
the Einstein relation. Second, a temperature gradient 



would not induce an intrinsic charge Hall current, violat- 
ing the Mott relation. Finally, it is also unclear whether 
the Onsager relation is satisfied or not. 

It turns out the correct description of anomalous ther- 
moelectric transport in ferromagnets requires the knowl- 
edge of both the magnetic moment and orbital magne- 
tization. Firstly, as we showed in Eq. (4.19), the local 
current is given by 



jL = / (Sr g(r > <^ + v x / jSy /(r < q)m{q) ' 

(4.21) 

where in the second term we have replaced the distri- 
bution function g(r, q) with the local Fermi-Dirac func- 
tion f(r,q), which is sufficient for a first-order calcula- 
tion. Secondly, in ferromagnetic systems, it is impor- 
tant to discount the contribution from the magnetiza- 
tion current. It was argued that the magnetization cur- 
rent cannot be measured by conventional transport ex- 
periments (Cooper et al., 1997). Therefore the transport 
current is given by 

j=j L -Vx M(r) . (4.22) 

Using Eq. (4.14), one finds 

3 = J ^f9{r,q)r-\vx J def(e)a^(e) . (4.23) 

Equation (4.23) is the most general expression for 
the transport current. We notice that the contribution 
from the orbital magnetic moment m(q) cancels out. 
This agrees with the intuitive picture we developed in 
Sec. IV. B, i.e., the orbital moment is due to the self- 
rotation of the wave packet, therefore it is localized and 
cannot contribute to transport (see Fig. 10). 

In the presence of a statistical force, there are two 
ways for a Hall current to occur. The asymmetric 
scattering will have an effect on the distribution g(r, q), 
which is obtained from the Boltzmann equation (Berger, 
1972). This results in a transverse current in the first 
term of Eq. (4.23). In addition to that, there is an 
intrinsic contribution comes from the orbital magneti- 
zation, which is the second term of Eq. (4.23). Note 
that the spatial dependence enters through T(r) and 
/i(r) in the distribution function. It is straightforward 
to verify that for the intrinsic contribution to the 
anomalous thermoelectric transport, both the Einstein 
relation and Mott relation still hold (Onoda et al., 2008; 
Xiao et al., 2006b). Hence, the measurement of this 
type of transport, such as the anomalous Nernst effect, 
can give further insight of the intrinsic mechanism of 
the anomalous Hall effect. Great experimental efforts 
have been put along this line. The intrinsic contribu- 
tion has been verified in CuC^Se^zBix (Lee et al., 
2004a,b), Lai_3.Sr3.CoO3 (Miyasato et al, 2007), 
Nd 2 Mo 2 7 and Sm 2 Mo 2 07 (Hanasaki et al, 2008), 
Gai-.Mn, (Pu et al, 2008). 

Equation (4.23) is not limited to transport driven by 
statistical forces. As we shall show later, at the micro- 
scopic level the mechanical force generally has two effects: 
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it can drive the electron motion directly, and appears in 
the expression for r; it can also make the electron energy 
and the Berry curvature spatially dependent, hence also 
manifest in the second term in Eq. (4.23). The latter 
provides another route for the Berry phase to enter the 
transport problems in inhomogcneous situations, which 
can be caused by a non-uniform distribution function, or 
a spatially-dependent perturbation, or both. 



V. ELECTRON DYNAMICS IN THE PRESENCE OF 
ELECTROMAGNETIC FIELDS 

In last section we discussed the construction and 
general properties of a wave packet. Now we are set to 
study its dynamics under external perturbations. The 
most common perturbations to a crystal is the electro- 
magnetic fields. The study of the electron dynamics 
under such perturbations dates back to Bloch, Peierls, 
Jones, and Zener in the early 1930s, and is continued 
by Slater (1949), Luttinger (1951), Adams (1952), 
Karplus and Luttinger (1954), Kohn and Luttinger 
(1957), Adams and Blount (1959), Blount (1962a), 
Brown (1967), Zak (1977), Rammal and Bcllissard 
(1990), Wilkinson and Kay (1996). In this section we 
present the semiclassical theory based on the wave 
packet approach (Chang and Niu, 1995, 1996). 



A. Equations of motion 

In the presence of electromagnetic fields, the Hamilto- 
nian is given by 

H= \p + eA(r)]* 
Zm 

where V(r) is the periodic lattice potential, and A(r) 
and 4>(r) are the electromagnetic potentials. If the length 
scale of the perturbations is much larger than the spatial 
spread of the wave packet, the approximate Hamiltonian 
that the wave packet "feels" may be obtained by lineariz- 
ing the perturbations about the wave packet center r c as 



AH 



H « H c + AH , 
[p+eA(r c )} 2 



(5.2) 

2w . .„ ™ (5.3) 

" {A(r) - A{r c ),p} - eE ■ (r - r c ) , (5.4) 



+ V(r) - ecj){r c ) , 



2m 



where {, } is the anticommutator. Naturally, we can then 
construct the wave packet using the eigenstates of the 
local Hamiltonian H c . The effect of a uniform A(r c ) 
is to add a phase to the eigenstates of the unperturbed 
Hamiltonian. Therefore the wave packet can be written 
as 



\W(k c ,r c )) = e-^ hAir ^ r \W (k c ,r c )) , 



(5.5) 



where |Wq) is the wave packet constructed using the un- 
perturbed Bloch functions. 



The wave packet dynamics can be obtained 
from the time-dependent variational princi- 
ples (Kramer and Saraceno, 1981). The basic recipe is to 
first obtain the Langragian from the following equation, 



H\W) 



(5.6) 



then obtain the equations of motion using the Eu- 
ler equations. Straightforward calculation shows that 
{W\ih§- t \W) = eA R c - h§t arg w{k c ,t). For the wave 
packet energy, we have (W|Ai?|H^) = —m(k) ■ B. This 
is expected as we already showed that the wave packet 
carries an orbital magnetic moment m(k) that will cou- 
ple to the magnetic field. Using Eq. (4.5), we find the 
Langragian is given by, up to some unimportant total 
time-derivative terms (dropping the subscript c on r c and 

fee), 

L = Tik ■ r - ejw(fe) + ecf)(r) - er ■ A(r,t) + Tik ■ An(k) , 

(5.7) 

where £A/(fe) = £o(k) — B ■ m(k) with £o(fe) being the 
unperturbed band energy. The equations of motion is 

r=^±-kxn(k), (5.8a) 

Tik = -eE -erxB. (5.8b) 

Compared to the conventional equations of motion for 
Bloch electrons (Ashcroft and Mermin, 1976), there are 
two differences: (1) The electron energy is modified by 
the orbital magnetic moment; (2) the electron velocity 
gains an extra velocity term proportional to the Berry 
curvature. As we can see, in the case of only an elec- 
tric field, Eq. (5.8a) reduces to the anomalous velocity 
formula (3.6) we derived before. 



B. Modified density of states 

The Berry curvature not only modifies the electron dy- 
namics, but also has a profound effect on the electron 
density of states in the phase space (Xiao et al. , 2005) . 

Recall that in solid state physics, the expectation value 
of an observable, in the Bloch representation, is given by 



22 fnh(j>nk\0\lpnk 



(5.9) 



irk 



where is the distribution function. In the semiclas- 
sical limit, the sum is converted to an integral in the 
fe-space, 



E 



dk 



(5.10) 



where V is the volume, and (27r) rf is the density of states, 
i.e., number of states per unit fc-volume. From a classical 
point of view, the constant density of states is guaranteed 
by the Liouvilles theorem, which states that the volume 
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element is a conserved quantity during the time evolu- 
tion of the system. 7 However, as we shall show below, 
this is no longer the case for the Berry-phase modified 
dynamics. 

The time evolution of a volume element AV = ArAk 
is given by 



1 dAV 
AV dt 



V r - r + V fc k 



(5.11) 



Insert the equations of motion (5.8) into the above equa- 
tion. After some algebra, we find 



AV = 



Vn 



(1 + (e/h)B ■ ft) 



(5.12) 



The fact that the Berry curvature is generally k- 
dependent and the magnetic fields is r-dependent im- 
plies that the phase-space volume AV" changes during 
time evolution of the state variables (r, k). 

Although the phase space volume is no longer con- 
served, it is a local function of the state variables and 
has nothing to do with the history of time evolution. We 
can thus introduce a modified density of states 



D(r,k) 



(2*) 



l - d {i+ e - B -n) 



(5.13) 



such that the number of states in the volume element, 
D n (r,k)AV, remains constant in time. Therefore, the 
correct semiclassical limit of the sum in Eq. (5.9) is 



0(R) = J dkD(r,k){OS(r-R)) w , 



(5.14) 



where (■ ■ -)w is the expectation value in a wave packet, 
which could includes the dipole contribution due to the 
finite size of the wave packet (See Sec. 4.18). In a uniform 
system it is simply given by 



= J dkD(k)f(k)0(k) 



(5.15) 



We emphasize that although the density of states is 
no longer a constant, the dynamics itself is still Hamil- 
tonian. The modification comes from the fact that the 
dynamical variables, r and fc, are no longer canonical 
variables, and the density of states can be regarded as 
the phase space measure (Bliokh, 2006b; Duval et ai, 
2006b, a; Xiao et ai, 2006a). A more profound reason 
for this modification has its quantum mechanical ori- 
gin in non-commutative quantum mechanics, discussed 
in Sec. VII. 

In the following we discuss two direct applications of 
the modified density of states in metals and in insulators. 



1. Fermi volume 

We show that the Fermi volume can be changed lin- 
early by a magnetic field when the Berry curvature is 
nonzero. Assume zero temperature, the electron density 
is given by 



dk 

(2tt) c 



1+-B fl )e(e F -e) 



(5.16) 



We work in the canonical ensemble by requiring the elec- 
tron number fixed, therefore, to first order of B, the 
Fermi energy must be changed by 



5V F = - 



dk e 
(2Tr) d h 



b n 



(5.17) 



It is particularly interesting to look at insulators, where 
the integration is limited to the Brillouin zone. Then the 
electron must populate a higher band if J BZ dk B ■ Q is 
negative. When this quantity is positive, holes must ap- 
pear at the top of the valance bands. Discontinuous be- 
havior of physical properties in a magnetic field is there- 
fore expected for band insulators with a nonzero integral 
of the Berry curvatures (Chern numbers) . 



2. Streda Formula 

In the context of the quantum Hall effect, Streda 
(1982) derived a formula relating the Hall conductivity 
to the field derivative of the electron density at a fixed 
chemical potential 



(5.18) 



There is a simple justification of this relation by a ther- 
modynamic argument by considering the following adia- 
batic process in two dimensions. A time dependent mag- 
netic flux generates an electric field with an emf around 
the boundary of some region; and the Hall current leads 
to a net flow of electrons across the boundary and thus 
a change of electron density inside. Note that this ar- 
gument is valid only for insulators because in metals the 
adiabaticity would break down. Using Eq. (5.16) for an 
insulator, we obtain, in 2D 



>xy 



dk 

h ,/bz W? aXV 



(5.19) 



This is what Thouless et al. (1982) obtained using the 
Kubo formula. The fact that the quantum Hall conduc- 
tivity can be derived using the modified density of states 
further confirms the necessity to introduce this concept. 



C. Orbital magnetization: Revisit 

The actual value of this constant volume for a quantum state, 

however, can be determined only from the quantization condi- We have discussed the Orbital magnetization using a 

tions in quantum mechanics. rather pictorial derivation in Sec. IV. B. Here we derive 
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the formula again by using the field-dependent density of 
states (5.13). 

The equilibrium magnetization density can be ob- 
tained from the grand canonical potential, which, within 
first order in the magnetic field, may be written as 

F = -i ^log(l + e-P^-ri) 
" k 



If dk , e „ „, , 



(5.20) 

where the electron energy em — e(k) — m(k) ■ B in- 
cludes a correction due to the orbital magnetic mo- 
ment m(k). The magnetization is then the field deriva- 
tive at fixed temperature and chemical potential, M = 
-(8F/dB)^ T , with the result 



/dk 



dk e 
(3 J (27r) d h 



n(fc)iog(i 



-/3(e-Mh 



(5.21) 

Integration by parts of the second term will give us the 
exact formula obtained in Eq. (4.14). We have thus 
derived a general expression for the equilibrium orbital 
magnetization density, valid at zero magnetic field but 
at arbitrary temperatures. From this derivation we can 
clearly see that the orbital magnetization is indeed a bulk 
property. The center-of-mass contribution identified be- 
fore comes from the Berry-phase correction to the elec- 
tron density of states. 

Following the discussions on band insulators in our first 
example in Sec. V.B.I, there will be a discontinuity of the 
orbital magnetization if the integral of the Berry curva- 
ture over the Brillouin zone, or the anomalous Hall con- 
ductivity, is non-zero and quantized. Depending on the 
direction of the field, the chemical potential fj,o in the 
above formula should be taken at the top of the valence 
bands or the bottom of the conduction bands. The size 
of the discontinuity is given by the quantized anomalous 
Hall conductivity times E g /e, where E g is the energy 
gap. 

Similar formula for insulators with zero Chcrn 
number has been obtained by Ceresoli et al. (2006); 
Thonhauscr et al. (2005) using the Wannier function ap- 
proach, and has been numerically confirmed for a tight- 
binding model. Recently, Shi et al. (2007) provided 
a full quantum mechanical derivation of the formula, 
and showed that it is valid in the presence of electron- 
electron interaction, provided the one-electron energies 
and wave functions are calculated sclf-consistently within 
the framework of the exact current and spin-density func- 
tional theory (Vignale and Rasolt, 1988). 

The appearance of the Hall conductivity is not a coin- 
cidence. Let us consider an insulator. The free energy is 
given by 



Using the Maxwell relation, we have 
/ dn\ ( dM \ 



B.T 



(5.23) 



On the hand, the zero-temperature formula of the mag- 
netization for an insulator is given by 

f dk p 
M =\ 7^ WW + rCM -£)«}. (5.24) 
Jbz (2tt) j h 

Inserting it into Eq. (5.23) gives us once again the quan- 
tized Hall conductivity. 



D. Magnetotransport 

The equations of motion (5.8) and the density of 
states (5.13) gives us a complete description of the elec- 
tron dynamics in the presence of electromagnetic fields. 
In this subsection we apply these results to the prob- 
lem of magnetotransport. For simple notations, we set 
e — Ti — 1 and introduce the shorthand [dk] = dk/(2n) d . 



1. Cyclotron period 

Scmiclassical motion of a Bloch electron in a uni- 
form magnetic field is important to understand various 
magneto-effects in solids. In this case, the equations of 
motion reduce to 

D(k)r = v + (v ■ fl)B , (5.25a) 
D(k)k = —v X B , (5.25b) 

where D(k) = D{k)/{2Ti) d = 1 + (e/K)B ■ ft. 

We assume the field is along the z-axis. From the sec- 
ond equation of (5.25) we can see that motion in fe-space 
is confined in the xy-planc and is completely determined 
once the energy e and the z component of the wave vector 
k z is given. Let us calculate the period of the cyclotron 
motion. The time for the wave vector to move from k\ 
to k 2 is 



dt 



fe 2 
fci 



dk 
Ifel 



From the equations of motion (5.25) we have 
... B\v ± \ B\(de/dk) ± \ 



D(k) 



TiD{k) 



(5.26) 



(5.27) 



On the other hand, the quantity (de/dk)± can be written 
as Ae/Ak, where AA; denotes the vector in the plane 
connecting points on neighboring orbits of energy e and 
e + Ae, respectively. Then 



dF = —MdB - nd/j, - SdT 



(5.22) 



t 2 - U = 



D{k)Akdk 
Ae 



(5.28) 
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Introducing the 2D electron density for given e and k z 



n 2 {e,k z ) 



D{k) dk x dk y 



lk z ,e(k)<e ( 27r ) 2 

the period of a cyclotron motion can be written as 



(5.29) 



T= (2tt 



2 h dn 2 (e, k z ) 



D 



de 



(5.30) 



We thus recovered the usual expression for the cyclotron 
period, with the 2D electron density, Eq. (5.29), defined 
with the modified density of states. 

In addition, we note that there is a velocity term pro- 
portional to B in Eq. (5.25), which seems to suggest there 
will be a current along the field direction. We show that 
after averaging over the distribution function, this cur- 
rent is actually zero. The current along B is given by 



3b = 



-eB J [dk]fv ■ ft 

-^bJ [dk]v k F-n 

-^B(J{dk]V k (Fn) 



[dk}FV k ■ n 



(5.31) 



where F(e) = - f(e')de' and f(e) = dF/de. The 
first term vanishes 8 and if there is no magnetic monopole 
in fc-space, the second term also vanishes. In above cal- 
culation we did not consider the change of the Fermi 
surface. Since it always comes in the form (df /d/i)Sfi = 
— {df/de)8^L we can use the same technique to prove that 
the corresponding current also vanishes. 



2. The high field limit 

We now consider the magnetotransport at the so-called 
high field limit, i.e., uj c t 3> 1, where u c — 2ir/T is the 
cyclotron frequency and r is the relaxation time. We 
consider configuration where the electric and magnetic 
fields are perpendicular to each other, i.e., E — Ex, B — 
Bz and E B = 0. 

In the high ficd limit, w c t» 1, the electron can finish 
several turns between two successive collisions. We can 
then assume all orbits are closed. According to the theo- 
rem of adiabatic drifting (Niu and Sundaram, 2001), an 



8 For any periodic function F(k) with the periodicity of a recip- 
rocal Bravais lattice, the following identity holds for integrals 
taken over a Brillouin zone, f BZ dk Vfc-F(fc) = 0. To see this, 
consider I(k') = J dkF(k + fc'). Because -F(fc) is periodic 
in fc, I(k') should not depend on k' . Therefore, Vfc//(fc') = 
J dk, V fe /F(fc + fc') = / dfcV fc F(fc + fc') = 0. Setting fc' = 
gives the desired expression. This is also true if F(k) is a vector 
function. 



originally closed orbit remains closed for weak perturba- 
tions, so that 







(k) = E + (r) x B 



Or 



E x B 

B 2 



(5.32) 



(5.33) 



The Hall current is simply the sum over (r)± of occupied 
states: 



3h 



Ex B 
E x B 



[dk]f(k)(l + B - ft) 
[dk]f(k)D(k) . 



(5.34) 



Therefore in the high field limit we reach the remarkable 
conclusion: the total current in crossed electric and mag- 
netic fields is the Hall current as if calculated from free 
electron model 



E x B 

-e 7: — n 

B 2 



(5.35) 



and it has no dependence on the relaxation time r. This 
result ensures that even in the presence of anomalous 
Hall effect, the high field Hall current gives the "real" 
electron density. 

Let us now consider the hole-like band. The Hall 
current is obtained by substracting the contribution of 
holes from that of the filled band, which is given by 
-eEx f[dk]Sl. Therefore 



jhole = e _ 



Ex B 



B 2 



[dk]D(k)[l - /(fc)] -eEx 



J [dk]U 



(5.36) 

So for the hole-like band, there is an additional term in 
the current expression proportional to the Chern number 
(the second integral) of the band. 



3. The Low Field Limit 

Next we consider the magnetotransport at the low field 
limit, i.e., oj c t <C 1. In particular, we show that the Berry 
phase induce a linear magnetoresistance. By solving the 
Boltzmann equation, one finds that the diagonal element 
of the conductivity is given by 



[dk\ 



9fo vj 
de D(k) 



(5.37) 



This is just the zeroth order expansion based on lo c t. 
There are four places in this expression depending on B. 
(1) There is an explicit B-dependence in D(k). (2) The 
electron velocity v x is modified by the orbital magnetic 
moment: 



1 d(e - m z B) 
ti dk x 



,(o) 



1 dm z 



B . 



(5.38) 
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(3) There is also a modification to the Fermi energy, given 
by Eq. (5.17). (4) The relaxation time r can also depend 
on B. In the presence of the Berry curvature, the collision 
term in the Boltzmann equation is given by 



df 

at 



coll 



= - / [dk']D(k')W kk , [f(k) - f(k')} , (5.39) 



where is the transition probability from k! to k 

state. In the relaxation-time approximation we make the 
assumption that a characteristic relaxation time exists so 
that 

f —^ = D(k) J[dk']B^W kk i[f(k)-f(k')}. (5.40) 

If we assume fi(fc) is smooth and Wkk 1 is localized, the 
relaxation time can be approximated by 



TQ 

D(k) 



r (l-|-B-n). (5.41) 



More generally, we can always expand the relaxation time 
to first order of (e/~h)B ■ f2, 



T = T + Tl-B ■ $7 , 

n 



(5.42) 



where t\ should be regarded as a fitting parameter within 
this theory. 

Expand expression (5.37) to first order of B, and take 
the spherical band approximation, we obtain 



a« = e 2 r B 

e 



,df (2ett z 2 2 



dm. 



% dk x x 



((m)"V / [dk]fn 



(5.43) 



where M is the effect mass tensor. The zero-field con- 
ductivity takes the usual form 



r(°) 



The ratio —(T^x/o'xx will then give us the linear magne- 
toresistance. 



-e^J[dk] d -^vl 



(5.44) 



A. Equations of motion 

We consider a slowly perturbed crystal whose Hamil- 
tonian can be expressed in the following form 



H[r,p;px(r,t),...p g (r,t)] , 



(6.1) 



where {/3i(r,t)} are the modulation function charac- 
terizing the perturbations. They may represent either 
gauge potentials of electromagnetic fields, atomic dis- 
placements, charge or spin density waves, helical mag- 
netic structures, or compositional gradients. Following 
the same procedure as we have done in last section, we 
expand the Hamiltonian around the wave packet center, 
and obtain 



H = H c + AH , 
H c = H[r,p; {A(r OJ t)}] , 



(6.2) 
(6.3) 



AH = J2 V re A(r c , t) ■ {(r - r c ), — } . (6.4) 



Since the local Hamiltonian H c maintains periodicity of 
the unperturbed crystal, its eigenstates take the Bloch 
form 

H c (r c , t)\ip q (r c , t)) = e c (r c , q, t)\ip q (r c , t)) , (6.5) 

where q is the Bloch wave vector and e c (r c ,q,t) is the 
band energy. Here we have dropped the band index n for 
simple notations. 

Following the discussion in Sec. I.D, we switch to the 
Bloch Hamiltonian H c (q,r c ,t) = e-'i r H c (r c ,t)e lqr , 
whose eigenstate is the periodic part of the Bloch func- 
tion, \u(q,r c ,t)) — e^ tq " r \il}{q,r Cl t)). The Berry vector 
potentials can be defined for each of the coordinates of 
the parameter space (q,r c ,t); for example, 



A t = (u\id t \u) 



(6.6) 



After constructing the wave packet using the local 
Bloch functions \ip q (r c ,i)), one can apply the time- 
dependent variational principle to find the Langragian 
governing the dynamics of the wave packet: 

L = -e + q c ■ r c + q c ■ Aq + r c ■ Ar + A t , (6.7) 



VI. ELECTRON DYNAMICS IN THE PRESENCE OF 
GENERAL PERTURBATIONS 



Note that the wave packet energy e = e c + Ae has a 
correction Ae from AiJ, 



In this section we present the general theory of elec- 
tron dynamics in slowly perturbed crystals (Panati et al. , 
2003; Shindou and Imura, 2005; Sundaram and Niu, 
1999). As expected, the Berry curvature enters into 
the equations of motion and modifies the density of 
states. The difference is that one needs to introduce the 
Berry curvature defined in the extended parameter space 
(r,q,t). Two physical applications are considered: elec- 
tron dynamics in deformed crystals, and adiabatic cur- 
rent induced by inhomogeneity. 



fin rin 
As = (W\AH\W) = -3?(ipi| • (e c - H c )\^-) . (6.8) 

or c oq 

From the Lagrangian (6.7) we obtain the following equa- 
tions of equation: 

ds < > «-> 

r c = (tt qr ■ r c + tt qq ■ q c ) - U qt , (6.9a) 

uqc 

de ~ «-> 
Qc = - + (tt rr ■ r c + tt xq ■ q c ) + U rt (6.9b) 
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where f2's are the Berry curvatures. For example, 



8 A 



d r/3 A ka 



(6.10) 



In the following we will also drop the subscript c on r c 
and q c . 

The form of the equations of motion is quite symmet- 
rical with respect to r and q, and there are Berry cur- 
vatures between every pair of phase space variables plus 
time. The term tt qt was identified as the adiabatic veloc- 
ity vector in Sec. II. In fact, if the perturbation is uniform 
in space (has the same period as the unperturbed crystal) 
and only varies in time, all the spatial derivatives vanish; 
we obtain 



de 
dq 



ft 



qi 







(6.11) 



The first equation is the velocity formula (2.5) obtained 

in Sec. II. The term Q qq was identified as the Hall con- 
ductivity tensor. In the presence of electromagnetic per- 
turbations, we have 



H = H [q + eA(r)]-ecf>{r > t) 



(6.12) 



Hence the local basis can be written as \u(r c ,q)) = 
\u(k)), where k = q + eA(r). One can verify that by us- 
ing the chain rule d Qa — dk a and d Ta — {d ra Ap)dk l} , Ae 
given in Eq. (6.8) becomes —m(k) ■ B, and the equations 
of motion (6.9) reduces to Eq. (5.8). The physics of quan- 
tum adiabatic transport and the quantum and anomalous 
Hall effect can be described from a unified point of view. 
The Berry curvature f2 rt plays a role like the electric 
force. The antisymmetric tensor £l rr is realized in terms 
of the magnetic field in the Lorenz force and is also seen 
in the singular form (<5-function like distribution) of dis- 
locations in a deformed crystal (Bird and Preston, 1988). 
Finally, the Berry curvature between r and q can be re- 
alized in deformed crystals as a quantity proportional to 
the strain and the electronic mass renormalization in the 
crystal (Sundaram and Niu, 1999). 



B. Modified density of states 

The electron density of states is also modified by the 
Berry curvature. Let us consider the time-independent 
case. To better appreciate the origin of this modification, 
we introduce the phase space coordinates £ = (r, q). The 
equations of motion can be written as 

IVe/3=V ?a £, (6.13) 

where T = fl — J is an antisymmetric matrix with 



n = 



1 f*f* 










n 1 ' 







(6.14) 



According to standard theory of Hamiltonian dynam- 
ics (Arnold, 1978), the density of states, which is pro- 
portional to the phase space measure is given by 



D{r,q) = 



1 



(2tt) c 



det(Q - J) 



(6.15) 



One can show that in the time-dependent case D(r,q) 
has the same form. 

Let us consider the following situations, (i) If the per- 
turbation is electromagnetic field, by the variable substi- 
tution k = q + eA(r), Eq. (6.15) reduces to Eq. (5.13). 
(ii) In many situations we arc aiming at a first-order cal- 
culation in the spatial gradient. In this case, the density 
of states given by 



D 



1 



(27T)- 



(6.16) 



Note that if the Berry curvature vanishes, Eq. (6.13) 
becomes the canonical equations of motion for Hamilto- 
nian dynamics, and r and q are called canonical vari- 
ables. The density of states is a constant in this case. 
The presence of the Berry curvature renders the vari- 
ables non-canonical and, as a consequence, modifies 
the density of states. The non-canonical variables is 
a common feature of Berry-phase participated dynam- 
ics (Littlejohn and Flynn, 1991). 

To demonstrate the modified density of states, 
we again consider the Rice-Mele model discussed in 
Sec. II.C.l. This time we introduce the spatial depen- 
dence by letting the dimerization parameter 5(x) vary in 
space. Using Eq. (1.19) we find 



At sin 2 \dj 



4(A 2 + i 2 cos 2 f+<5 2 sin 2 !) 3 / 2 



(6.17) 



At half filling, the system is an insulator and its electron 
density is given by 



n dq 
.,2tt 



(6.18) 



We let 5(x) have a kink in its profile. Such a domain wall 
is known to carry fractional charge (Rice and Mele, 1982; 
Su et ai, 1979). Figure 11 shows the calculated electron 
density using Eq. (6.18) together with numerical result 
obtained by direct diagonalization of the tight-binding 
Hamiltonian. These two results are virtually indistin- 
guishable in the plot, which confirms the Berry-phase 
modification to the density of states. 



C. Deformed Crystal 

In this subsection we present a general the- 
ory of electron dynamics in crystals with deforma- 
tion (Sundaram and Niu, 1999), which could be caused 
by external pressure, defects in the lattice, or interfacial 
strain. 



29 



1.014 h 
1.012 
1.01 

'3 1.008 

G 

<d 

Q 1.006 

1.004 [■ 
1.002 
1 



■ Numerics 



Theory 




100 200 300 400 500 



J V. 



100 



200 300 

Lattice 



400 



500 



FIG. 11 (color online) Electron density of the Rice-Mele 
model with a spatial varying dimcrization parameter. The 
parameters we used are A = 0.5, t = 2, and 8 = tanh(0.02:r). 
Inset: The profile of S(x). From Xiao et al. 



Let us set up the basic notations for this problem. 
Consider a deformation described by the atomic displace- 
ment {ui}. We denote the deformed crystal potential as 
V(r; {R[ + Ui}), where R[ is the atomic position with I 
labeling the atomic site. Introducing a smooth displace- 
ment field u(r) such that u(Ri +Ui) — it;, the Hamilto- 
nian can be written as 



H 



P 

2m 



V[r - u(r)] + s aP (r)V a0 [r - u(r)] , (6.19) 



where s a p — du a /drp is the unsymmetrized strain, and 
V a p[r - u(r)] = ^[Rt + u(r) - r]p(dV/dRi a ) is a gra- 
dient expansion of the crystal potential. The last term, 
being proportional to the strain, can be treated pertur- 
batively. The local Hamiltonian is given by 



H c = f- + V[r - u(r c )} , 
Ira 



(6.20) 



with its eigenstates \ip q (r — u(r c ))). 

To write down the equations of motion, two pieces of 
information are needed. One is the gradient correction 
to the electron energy, given in Eq. (6.8). It is found 
that (Sundaram and Niu, 1999) 



where 



Ae = s a/3 D a/3 (q) 



D a/ 3 = m[v a Vp - (v a Vp)] + (V a p) 



(6.21) 



(6.22) 



with (• • •) the expectation value of the enclosed operators 
in the Bloch state, and v a is the velocity operator. Note 
that in the free electron limit (V — > 0) this quantity 
vanishes. Which is expected since a wave packet should 
not feel the effect of a deformation of the lattice in the 
absence of electron-phonon coupling. The other is the 
Berry curvature, which is derived from the Berry vector 



potentials. For deformed crystals, in addition to A. q , 
there are two other vector potentials 



ft 



with 



du a 
dr 



f(q) 



fa 



du a 

Ik 



to de 
h dq 



fiq 



(6.23) 



(6.24) 



It then leads to the following Berry curvatures 

du-f df-y du-y 9/ 7 



drp dq a 



Qk a t — 



dt dq a ' (6.25) 



With the above information we just need plug in the elec- 
tron energy as well as the Berry curvatures into Eq. (6.9) 
to obtain the equations of motion. 

We first consider a one-dimensional insulator with lat- 
tice constant a and is under a uniform strain with a new 
lattice constant a + da, i.e., d x u — Sa/a. Assuming one 
electron per unit cell, the electron density goes from 1 /a 
to 



1 



da a 



da 



)• 



(6.26) 



On the other hand, we can also directly calculate the 
change of the electron density using the modified density 
of states (6.16), which gives 



2tt / a j c 

dq da 



2tt 



(6.27) 



From a physical point of view, it says an insulator under 
a uniform strain remains an insulator. 

The above formalism is also applicable to dislocation 
strain fields, which are well defined except in a region 
of a few atomic spacings around the lin of dislocation. 
Outside this region, the displacement field u(r) is a 
smooth but multiple-valued function. On account of this 
multiple- valucdness, a wave packet of incident wave vec- 
tor q taken around the line of dislocation acquires a Berry 
phase 



7= (j)dr-Ar = jdu- f(k) « b ■ f(k) 



(6.28) 



where b — §dr a du/dr a is known as the Burgers vec- 
tor. What we have here is a situation similar to 
the Aharonov-Bohm effect (Aharonov and Bohm, 1959), 
with the dislocation playing the role of the solenoid, and 
the Berry curvature tt rr the role of the magnetic field. 
Bird and Preston (1988) showed that this Berry phase 
can affect the electron diffraction pattern of a deformed 
crystal. 

The above discussion only touches a few general ideas 
of the Berry phase effect in deformed crystals. With 
the complete information of the equations of motion, the 
scmiclassical theory provides a powerful tool to investi- 
gate the effects of deformation on electron dynamics and 
equilibrium properties. 
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D. Polarization induced by inhomogeneity 

In Sec. II. C we have discussed the Berry phase theory 
of polarization in crystalline solids, based on the basic 
idea that the polarization is identical to the integration 
of the adiabatic current flow in the bulk. There the sys- 
tem is assumed to be periodic and the perturbation de- 
pends only on time (or any scalar for that matter). In 
this case, it is straightforward to obtain the polarization 
based on the equations of motion (6.11). However, in 
many physical situations the system is in an inhomogc- 
neous state and the electric polarization strongly depends 
on the inhomogeneity. Examples include flexoelectricity 
where a finite polarization is produced by a strain gra- 
dient (Tagantscv, 1986, 1991), and multiferroic materials 
where the magnetic ordering varies in space and induces a 
polarization (Cheong and Mostovoy, 2007; Fiebig et al, 
2002; Hur et al, 2004; Kimura et al, 2003). 

Let us consider an insulating crystal with an order pa- 
rameter that varies slowly in space. We assume that, 
at least at the mean-field level, the system can be de- 
scribed by a perfect crystal under the influence of an 
external field h(r). If, for example, the order parame- 
ter is the magnetization, then h(r) can be chosen as the 
exchange field that yields the corresponding spin config- 
uration. Our goal is to calculate the electric polarization 
to first order in the spatial gradient as the field h(r) is 
gradually turned on. The Hamiltonian thus takes the 
form H[h(r); A] where A is the parameter describing the 
turning on process. Xiao et al. (2007a) showed that the 
first order contribution to the polarization can be clas- 
sified into two categories: the perturbative contribution 
due to the correction to the wave function, and the topo- 
logical contribution which is from the dynamics of the 
electrons. 

Let us first consider the perturbation contribution, 
which is basically a correction to the polarization formula 
obtained by King-Smith and Vandcrbilt (1993) for a uni- 
form system. The perturbative contribution is obtained 
by evaluating the Berry curvature U, qt in Eq. (2.27) to 
first order of the gradient. Remember that we always ex- 
pand the Hamiltonian into the form H = H c + AH, and 
choose the eigenfunctions of H c as our expansion basis. 
Hence to calculate the Berry curvature to first order of 
the gradient, one needs to know the form of the wave 
function perturbed by AH. It has been extensively dis- 
cussed in the case of an electric field (Nunes and Gonze, 
2001; Souza et al, 2002). 

The topological contribution is of different nature. 
Starting from Eq. (6.9) and making use of the modified 
density of states (6.16), one finds the adiabatic current 
induced by inhomogeneity is give by 

(2) _ f d( l ( Qqq QrX , Qqr Qq\ _ Qqr ^qX\ 

3 a ~ e J BZ (2-K) d \ aP (3 W a Q/3 13 )' 

(6.29) 

We can see that this current is explicitly proportional 
to the spatial gradient. Comparing this equation with 
Eq. (2.6) reveals a very elegant structure: the zeroth 



order contribution, Eq. (2.6), is given as an integral 
of the first Chern form, while the first order contribu- 
tion, Eq. (6.29), is given as an integral of the second 
Chern form. Similar result has been obtained by Qi et al. 
(2008). 

The polarization is obtained by integrating the current. 
As usually in the case of multiferroics, we can assume 
the order parameter is periodic in space (but in general 
incommensurate with the crystal lattice). A two-point 
formula can be written down 9 

(6.30) 

where V is the volume of the periodic structure of the 
order parameter. Again, due to the loss of tracking of 
A, there is an uncertain quantum which is the second 
Chern number. If we assume the order parameter has 
period l y in the y-direction, the polarization quantum in 
the x-direction is given by 



where a is the lattice constant. 

Kunz (1986) has discussed the charge pumping in in- 
commensurate potentials and he showed that in general 
the charge transport is quantized and is given in the form 
of Chern numbers, which is consistent with what we have 
derived. 

The second Chern form demands that the system 
must be two-dimensional or higher, otherwise the sec- 
ond Chern form vanishes. It allows us to determine the 
general form of the induced polarization. Consider a two- 
dimensional minimal model with h(r) having two com- 
ponents. If we write H[h(r);X] as H[\h(r)], i.e., A acts 
like a switch, the polarization can be written as 

PW = ^f dr X [(V ■ h)h - (h ■ V)h] . (6.32) 

The coefficient x 1S given by 

[ tttx^ f d\e abcd n ab n cd , (6.33) 
8 Jbz [ 2 V Jo 

where the Berry curvature is defined on the parameter 
space (q,h), and t a bcd is the Levi-Civita antisymmetric 
tensor. 

Xiao et al. (2007a) also showed how the two-point for- 
mula can be implemented in numerical calculations using 
a discrctized version (Kotiuga, 1989). 



9 So far we only considered the Abelian Berry case. The non- 
Abelian reslut is obtained by replacing the Chern-Simons form 
with its non- Abelian form 
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1. Magnetic field induced polarization 

Using the above formalism, Essin et al. (2008) consid- 
ered the polarization induced by a magnetic field, which 
can be regarded as arising from a spatial inhomogene- 
ity in the vector potential A(r). Consider A = Byz 
with B = h/ea z l y x. Note that magnetic flux over the 
supercell a z x l y in the x-direction is exactly h/e, there- 
fore the system is periodic in the y-direction with period 
l y . According to our discussion in Sec. V, the effect of a 
magnetic field can be counted by the Peierls substitution, 
k z — * k z + eBy/h, hence V y — (eB /K)V kz- The induced 
polarization is given by 

9e 2 , s 

= 2nh^ ' ^ 

with 

If 2 

9 = — / dke a 3~ Tr[A a daAy — i-AaApAy] . (6.35) 
2lT Jbz a 

Recall that P = d£/dE, such a magnetic-field-induced 
polarization implies that there is an electromagnetic cou- 
pling of the form 

A£ EM - ' B ■ ( 6 - 36 ) 

Ann 

This coupling, labeled "axion electrodynamics" , was dis- 
cussed by Wilczek (1987). When 9 — ir, the correspond- 
ing insulator is known as a 3D Zi topological insual- 
tor (Qi et ai, 2008). 

E. Spin Texture 

So far our discussion has focused on the physical effects 
of the Berry curvature in the momentum space (ilfefc), 
or in the mixed space of the momentum coordinates and 
some other parameters (£lkr and fifct). In this subsection 
we discuss the Berry curvatures which originate only from 
the nontrivial real space configuration of the system. 

One of such systems is magnetic materials with domain 
walls or spin textures. Let us consider a ferromagnetic 
thin films described by the following Hamiltonian 

P 2 

H = — - Jn(r,t) ■ cr , (6.37) 
2m 

where the first term is the bare Hamiltonian for a con- 
duction electron, and the second term is the s-d cou- 
pling between the conduction electron and the local d- 
electron spin along the direction n(r,t) with J being 
the coupling strength. Note that we have allowed the 
spin texture to vary in both space and time. The simple 
momentum-dependence of the Hamiltonian dictates that 
all fc-dependent Berry curvatures vanish. 

Because of the strong s-d coupling, we adopt the adi- 
abatic approximation which states that the electron spin 
will follow the local spin direction during its motion. 



Then the spatial variation of local spin textures gives 
rise to the Berry curvature field 

tt rr = l - sin 9{V9 x Vcj)) , (6.38) 

where 9 and 4> are the spherical angles specifying the 
direction of h. According to Eqs. (6.9), this field acts on 
the electrons as an effective magnetic field. In addition, 
the time-dependence of the spin texture also gives rises 
to 

Sl rt = i sin 9{d t <j)V9 - d t 9V<j>) . (6.39) 

Similarly, Sl rt acts on the electrons as an effective electric 
field. This is in analogy with a moving magnetic field 
(r2 rr ) generating an electric field {Q. rt ). 

The physical consequences of these two fields are ob- 
vious by analogy with the electromagnetic fields. The 
Berry curvature tt rr will drive a Hall current, just like 
the ordinary Hall effect (Bruno et al., 2004; Ye et ai, 
1999). Unlike the anomalous Hall effect discussed in 
Sec. III.D, this mechanism for a nonvanishing Hall ef- 
fect does not require the spin-orbit coupling, but does 
need a topologically nontrivial spin texture, for example, 
a skyrmion. On the other hand, for a moving domain wall 
in a thin magnetic wire, the Berry curvature n rt will in- 
duce an electromotive force, which results in a voltage 
difference between the two ends. This Berry curvature 
induced emf has has been experimentally measured re- 
cently (Yang et ai, 2009). 



VII. QUANTIZATION OF ELECTRON DYNAMICS 

In previous sections, we have reviewed several Berry 
phase effects in solid state systems. Berry curvature of- 
ten appears as a result of restricting (or projecting) the 
extent of a theory to its subspace. In particular, the 
Berry curvature plays a crucial role in the semiclassi- 
cal dynamics of electrons, which is valid under the one- 
band approximation. In the following, we will explain 
how could the semiclassical formulation be re-quantized. 
This is necessary, for example, in studying the quan- 
tized Wannier-Stark ladders from the Bloch oscillation, 
or the quantized Landau levels from the cyclotron orbit 
(Ashcroft and Mermin, 1976). The re-quantized theory 
is valid in the same subspace of the semiclassical theory. 
It will become clear that, the knowledge of the Bloch 
energy, the Berry curvature, and the magnetic moment 
in the semiclassical theory constitute sufficient informa- 
tion for building the re-quantized theory. In this section, 
we focus on the following methods of quantization: the 
Bohr-Sommerfeld quantization and the canonical quan- 
tization. 



32 



A. Bohr-Sommerfeld quantization 



B. Wannier-Stark ladder 



A method of quantization is a way to select quantum 
mechanically allowed states out of a continuum of clas- 
sical states. This is often formulated using the general- 
ized coordinates qt and their conjugate momenta pi. The 
Bohr-Sommerfeld quantization requires the action inte- 
gral for each set of the conjugate variables to satisfy 



Pidqi 



(m,i + ~J h. / 



1, 



(7.1) 



where Ci are closed loops in the phase space with dimen- 
sion 2d, rrii are integers, and Vi are the so-called Maslov 
indices, which are usually integers. 

However, since the choices of conjugate variables are 
not unique, the Bohr-Sommerfeld quantization method 
may give inconsistent quantization rules. This is known 
to happen in the case of an isotropic three-dimensional 
harmonic oscillator (Tabor, 1989). This problem can be 
fixed if, instead of Eq. (7.1), one uses the following inte- 
grals, 



« d 



h.k = 1, 



,d, (7.2) 



where Ck are the periodic orbits on invariant tori. An in- 
variant torus is a torus with constant action in the phase 
space. The revised rule above is often called the Einstein- 
Brillouin- Keller (EBK) quantization (Tabor, 1989). 

In the wavepacket formulation of Bloch electrons, 
both r c and q c are treated as generalized coordinates. 
With the Lagrangian in Eq. (5.7), one can find then- 
conjugate momenta dL/dr c and dL/dq c , which are 
equal to hq c and Ti(u\idu/dq c ) — fiA respectively 
(Sundaram and Niu, 1999). The quantization condition 
for an orbit with constant energy thus becomes, 



2vr ra + 



2tt 



(7.3) 



where Tc = § C A ■ dq c is the Berry phase of a constant- 
energy orbit C (also see Kuratsuji and Iida (1985); 
Wilkinson (1984b)). Since the Berry phase is path depen- 
dent, one may need to solve the equation self-consistently 
to obtain the quantized orbits. 

Before applying the Bohr-Sommerfeld quantization in 
the following subsections, we would like to point out 
two disadvantages of this method. First, the value of 
the Maslov index is not always evident. For exam- 
ple, for an one-dimensional particle bounded by two 
walls, its value would depend on the slopes of the walls 
(van Houten et ai, 1989). In fact, a non-integer value 
may give a more accurate prediction of the energy levels 
(Friedrich and Trost, 1996). Second, this method fails if 
the trajectory in phase space is not closed, or if the dy- 
namic system is chaotic so that no invariant torus exists. 
On the other hand, the method of canonical quantization 
in Sec. VII. D does not have these problems. 



Consider an electron moving in a one-dimensional pe- 
riodic lattice with lattice constant a. Under a weak 
uniform electric field E, according to the semiclassical 
equations of motion, the quasi-momentum of an electron 
wavepacket is simply (see Eq. (5.8)) 



hq c (t) = -eEt. 



(7.4) 



It takes the time Tg = h/(eEa) for the electron to tra- 
verse the first Brillouin zone. Therefore, the angular fre- 
quency of the periodic motion is lob = eEa/Ti. This is the 
so-called Bloch oscillation (Ashcroft and Mermin, 1976). 

Similar to a simple harmonic oscillator, the energy 
of the oscillatory electron is quantized in multiples of 
Tilob- However, unlike the former, the Bloch oscillator 
has no zero-point energy (that is, the Maslov index is 
zero). These equally spaced energy levels are called the 
Wannier-Stark ladders. Since the Brillouin zone is peri- 
odic, the electron orbit is closed. According to the Bohr- 
Sommerfeld quantization, one has 



dq c 



-27r m 



r. 



2tt 



(7.5) 



For a simple one-dimensional lattice with inversion sym- 
metry, the Berry phase Tc m can only have two values, 
or 7r (Zak, 1989b), as discussed earlier in Sec. II. C. 

Starting from Eq. (7.5), it is not difficult to find the 
average position of the electron, 



2tt 



(7.6) 



where we have neglected the subscript m in Tc m 
since all of the paths in the same energy band have 
the same Berry phase here. Such average positions 
{fcjm are the average positions of the Wannier function 
(Vanderbilt and King-Smith, 1993). Due to the Berry 
phase, they are displaced from the positive ions located 
at am. 

In Sec. II. C, the electric polarization is derived using 
the theory of adiabatic transport. It can also be ob- 
tained from the expectation value of the position op- 
erator directly. Because of the charge separation men- 
tioned above, the one-dimensional crystal has a polar- 
ization AP — eT c /2ir (compared to the state without 
charge separation), which is the electric dipole per unit 
cell. This is consistent with the result in Eq. (2.28). 

After time average, the quantized energies of the elec- 
tron wavepacket are, 



(Or 



(e(q c )) - eE(r c ) 



eo — eEa m 



cj m 

r 



c 
2tt 



(7.7) 



which are the energy levels of the Wannier-Stark ladders. 

Two short comments are in order: First, beyond the 
one-band approximation, there exist Zener tunnellings 
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between Bloch bands. Therefore, the quantized levels are 
not stationary states of the system. They should be un- 
derstood as resonances with finite life-times (Gliick et ai, 
1999). Second, the fascinating phenomenon of Bloch 
oscillation is not commonly observed in laboratory for 
the following reason: In an usual solid, the electron 
scattering time is shorter than the period Tg by sev- 
eral orders of magnitude. Therefore, the phase coher- 
ence of the electron is destroyed within a tiny fraction 
of a period. Nonetheless, with the help of a superlat- 
tice that has a much larger lattice constant, the period 
Tb can be reduced by two orders of magnitude, which 
could make the Bloch oscillation and the accompanying 
Wannier-Stark ladders detectable (Mendes and Bastard, 
1993). Alternatively, the Bloch oscillation and Wannier- 
Stark ladders can also be realized in an optical lattice 
(Ben Dahan et ai, 1996; Wilkinson and Kay, 1996), in 
which the atom can be coherent over a long period of 
time. 



C. de Haas-van Alphen oscillation 

When an uniform B field is applied to a solid, the 
electron would execute a cyclotron motion in both the 
r-space and the fc-space. From Eq. (5.25b), it is not dif- 
ficult to see that an orbit C in fc-space lies on the in- 
tersection of a plane perpendicular to the magnetic field 
and the constant-energy surface (Ashcroft and Mermin, 
1976). Without quantization, the size of an orbit is deter- 
mined by the initial energy of the electron and can be var- 
ied continuously. One then applies the Bohr-Sommerfeld 
quantization rule, as Onsager did, to quantize the size 
of the orbit (Onsager, 1952). That is, only certain or- 
bits satisfying the quantization rule arc allowed. Each 
orbit corresponds to an energy level of the electron (i.e., 
the Landau level). Such a method remains valid in the 
presence of the Berry phase. 

With the help of the semiclassical equation (see 
Eq. (5.8)), 

Kk c = -er c x B, (7.8) 

the Bohr-Sommerfeld condition in Eq. (7.3) can be writ- 
ten as (notice that ?iq c = Tik c — eA, and v = 2), 

*J c JcX d *c={ m +l-^t)^ (™) 

where 4>q = h/e is the flux quantum. The integral on the 
left-hand side is simply the magnetic flux enclosed by the 
real-space orbit (allowing a drift along the B-direction) . 
Therefore, the enclosed flux has to jump in steps of the 
flux quantum (plus a Berry phase correction) . 

Similar to the Bohr atom model, in which the elec- 
tron has to form a standing wave, here the total phase 
acquired by the electron after one circular motion also 
has to be integer multiples of 2ir. Three types of phases 
contribute to the total phase: (a), The Aharonov-Bohm 



phase: an electron circulating a flux quantum picks up a 
phase of 2ir. (b), The phase lag of it at each turning point 
(there are two of them). This explains why the Maslov 
index is two. (c), The Berry phase intrinsic to the solid. 
Therefore, Eq. (7.9) simply says that the summation of 
these three phases should be equal to 2itm. 

The orbit in fc-space can be obtained by re-scaling the 
r-space orbit in Eq. (7.9) with a linear factor of X B , 
followed by a rotation of ninety degrees, where Xb = 
y^h/eB is the magnetic length (Ashcroft and Mermin, 
1976). Therefore, one has 

|.£k c x dkc = 2.( m+ i- r g.)f . (7.10) 

The size of the orbit, combined with the knowledge of 
the electron energy -E(k) c = e(k c ) — M • B, help deter- 
mining the quantized energy levels. For an electron with 
a quadratic energy dispersion (before applying the mag- 
netic field), these levels are equally spaced. However, 
with the Berry phase correction, which are usually dif- 
ferent for different orbits, the energy levels are no longer 
uniformly distributed. This is related to the discussion 
in Sec. V on the relation between the density of states 
and the Berry curvature (Xiao et al, 2005). 

As a demonstration, we apply the quantization rule 
to graphene and calculate the energies of Landau lev- 
els near the Dirac point. Before applying a magnetic 
field, the energy dispersion near the Dirac point is linear, 
i?(k) = hvpk. It is known that, if the energy disper- 
sion near a degenerate point is linear, then the cyclotron 
orbit will acquire a Berry phase Tc = 7T, independent 
of the shape of the orbit (Blount, 1962b). As a result, 
the 1/2 on the right hand side of Eq. (7.10) is cancelled 
by the Berry phase term. According to Eq. (7.10), the 
area of a cyclotron orbit is thus nk 2 — 2n\m\^- ,m € Z. 
From which one can easily obtain the Landau level en- 
ergy E m — vp\/2eBhm. The experimental observation 
of a quantum Hall plateau at zero energy is thus a direct 
consequence of the Berry phase (Novoselov et al, 2005, 
2006; Zhang et al, 2005). 

In addition to point degeneracy, other types of degener- 
acy in momentum space can also be a source of the Berry 
phase. For example, the effect of the Berry phase gener- 
ated by a line of band contact on magneto-oscillations is 
studied in Mikitik and Sharlai (1999, 2004). 

The discussion so far is based on the one-band approx- 
imation. In reality, the orbit in one band would couple 
with the orbits in other bands. As a result, the Landau 
levels are broadened into mini-bands (Wilkinson, 1984a). 
Similar situation occurs in a magnetic Bloch band, which 
is the subject of Sec. VIII. 



D. Canonical quantization (Abelian case) 

In addition to the Bohr-Sommerfeld quantization, an 
alternative way to quantize a classical theory is by finding 
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out position and momentum variables that satisfy the 
following Poisson brackets, 

{Xi,pj} = 8ij. (7.11) 

Afterwards, these classical canonical variables are pro- 
moted to operators that satisfy the commutation rela- 
tion, 

[xi,pj]=i%5ij, (7-12) 

That is, all we need to do is to substitute the Poisson 
bracket {xi,pj} by the commutator [xi,pj]/ih. Based on 
the commutation relation, these variables can be written 
explicitly using either the differential-operator represen- 
tation or the matrix representation. Once this is done, 
one can proceed to obtain the eigenvalues and eigenstates 
of the Hamiltonian iJ(x,p). 

Even though one can always have canonical pairs in a 
Hamiltonian system, as guaranteed by the Darboux the- 
orem (Arnold, 1989), in practice, however, finding them 
may not be a trivial task. For example, the center-of- 
mass variables r c and k c in the semiclassical dynamics 
in Eq. (5.8) are not canonical variables since their Pois- 
son brackets are not of the canonical form (Duval et al. , 
2006a; Xiao et al, 2005), 

{n,rj} = eijkClk/K, (7-13) 
{h,kj} = -eij k eB k /K, (7-14) 
{r^kj} = (Sij + eBiflj) / k, (7.15) 

where «(k) = 1 + eB ■ f2(k). In order to carry out 
the canonical quantization, canonical variables of posi- 
tion and momentum must be found. 

Let us start with two special cases. The first is a solid 
with zero Berry curvature that is under the influence of a 
magnetic field (ft — 0, B ^ 0). In this case, the factor k 
in Eq. (7.14) reduces to one and the position variables 
commute with each other. Obviously, if one assumes 
Tik c = p + eA(x) and requires x and p to be canon- 
ical conjugate variables, then the quantized version of 
Eq. (7.14) (with ih inserted) can easily be satisfied. This 
is the familiar Peierls substitution (Pcierls, 1933). 

In the second case, consider a system with Berry cur- 
vature but not in a magnetic field (fJ ^ 0, B = 0). In this 
case, again we have k = 1. Now the roles of r c and k c 
in the commutators are reversed. The momentum vari- 
ables commute with each other, but not the coordinates. 
One can apply a Peierls-like substitution to the coordi- 
nate variables and write r c = x + A(q) . It is not diffi- 
cult to see that the commutation relations arising from 
Eq. (7.13) can indeed be satisfied. After the canonical 
quantization, x becomes id/d<\ in the quasi-momentum 
representation. In Blount (1962b), the position opera- 
tor r in the one-band approximation acquires a correc- 
tion, which is our Berry connection A. Therefore, r c can 
be identified with the projected position operator PrP, 
where P projects to the energy band of interest. 



When both B and Q are nonzero, applying both of the 
Peierls substitutions simultaneously is not enough to pro- 
duce the correct commutation relations, mainly because 
of the non-trivial factor k there. In general, exact canon- 
ical variables cannot be found easily. However, since the 
semiclassical theory itself is valid to linear order of field, 
we only need to find the canonical variables correct to the 
same order in practice. The result is (Chang and Niu, 
2008), 

r c = x + A(ir) + G(tt), 
Kk c = p + eA(x) + eB x A(ir), (7.16) 

where n = p + eA(x), and G Q (k c ) = (e/K)(A x B) • 
dA/dk ca . This is the generalized Peierls substitution for 
systems with Berry connection A and vector potential A. 
With these equations, one can verify Eqs. (7. 13), (7. 14), 
and (7.15) to linear orders of B and f2. 

A few comments are in order: First, if a physical ob- 
servable is a product of several canonical variables, the 
order of the product may become a problem after the 
quantization since the variables may not commute with 
each other. To preserve the hermitian property of the 
physical observable, the operator product needs to be 
symmetrized. Second, the Bloch energy, Berry curvature, 
and orbital moment of the semiclassical theory contains 
sufficient information for building a quantum theory that 
accounts for all physical effects to first order in external 
fields. We will come back to this in Sec. IX, where the 
non-Abelian generalization of the canonical quantization 
method is addressed. 



VIII. MAGNETIC BLOCH BAND 

The semiclassical dynamics in previous sections is valid 
when the external field is weak, such that the latter can 
be treated as a perturbation to the Bloch states. Such 
a premise is no longer valid if the external field is so 
strong that the structure of the Bloch bands is signifi- 
cantly altered. This happens, for example, in quantum 
Hall systems where the magnetic field is of the order of 
Tesla and a Bloch band would break into many subbands. 
The translational symmetry and the topological property 
of the subband are very different from those of the usual 
Bloch band. To distinguish between the two, the former 
is called the magnetic Bloch band (MBB). 

The MBB usually carries non-zero quantum Hall con- 
ductance and has a nontrivial topology. Compared to the 
usual Bloch band, the MBB is a more interesting play- 
ground for many physics phenomena. In fact, the Berry 
curvature of the Bloch electron is first revealed in the 
MBB. In this section, we review some basic fact of the 
MBB, as well as the semiclassical dynamics of the mag- 
netic Bloch electron when it is subject to further elec- 
tromagnetic perturbation. Such an elegant formulation 
provides a clear picture of the hierarchical subbands split 
by the strong magnetic field (the Hofstadter spectrum) 
(Hofstadter, 1976). 
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A. Magnetic translational symmetry 

In the presence of a strong magnetic field, one needs to 
treat the magnetic field and the lattice potential on equal 
footing and solve the following Schrodinger equation, 



plaquette 



Magnetic unit cell 



[p + eA(r)] 2 + Vl(t) \ ^(r) = £ty(r), (8.1) 



1 
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where Vl is the periodic lattice potential. For conve- 
nience of discussion, we assume the magnetic field is uni- 
form along the z-axis and the electron is confined to the 
x-y plane. Because of the vector potential, the Hamilto- 
nian H above no longer has the lattice translation sym- 
metry. 

However, since the lattice symmetry of the charge 
density is not broken by an uniform magnetic field, 
one should be able to define translation operators 
that differ from the usual ones only by phase factors 
(Lifshitz and Landau, 1980). First, consider a system 
translated by a lattice vector a, 

~ [p + eA(r + a)] 2 + F L (r)| ^(r + a) = E^(r + a), 

(8.2) 

where Vt(r + a) = Vt(r) has been used. One can write 

A(r + a) = A(r) + V/(r), (8.3) 

where V/(r) = A(r + a) — A(r) = AA(a). For an uni- 
form magnetic field, the vector potential must be a lin- 
ear function of r. Therefore, A A is independent of r 
and / = A A • r. The extra vector potential V/ can be 
removed by a gauge transformation, 
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p + eA(r)] 2 + y L (r) 



ip(r+a) = £e's j ^(r+a) 
(8.4) 

We now identify the state above as the magnetic trans- 
lated state T a i/>(r), 
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FIG. 12 (color online) When the magnetic flux per plaquette 
is <j)o /3, the magnetic unit cell is composed of three plaquettes. 
The magnetic Brillouin zone is three times smaller than the 
usual Brillouin zone. Furthermore, the magnetic Bloch states 
are three-fold degenerate. 



eigenstate of H , T qai , and T a2 is called a magnetic Bloch 
state, and its energy the magnetic Bloch energy, 



q&l Ip' 



T 9ai V„k = e JK ' 9ai (k)Vnk, 
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(8.7) 
(8.8) 
(8.9) 



Since the magnetic unit cell is q times larger than the 
usual unit cell, the magnetic Brillouin zone (MBZ) has to 
be q times smaller. If bi and b2 are defined as the lattice 
vectors reciprocal to ai and a 2 . Then, in this example, 
the MBZ is folded back q times along the bi direction. 

In addition, with the help of Eqs. (8.6) and (8.8), one 
can show that the eigenvalues of the T &2 operator for the 
following translated states, 



(8.10) 



T a V>(r) = e*K* A - r ^(r + a) 



(8.5) 



an; 



The operator T a being defined this way has the desired 
property that [H,T a ]=0. 

Unlike usual translation operators, magnetic transla- 
tions along different directions usually do not commute. 
For example, let ai and a 2 be lattice vectors, then 



t(k+b 2 p/t/)-a 2 i(k+2b 2 p/<j)-a 2 



= i(k+(<3-l)b 2 p/tj)-a2 



T &2 T ai =T ai T &2 e^ 



i A dr 



(8.6) 



where § A • dr is the magnetic flux going through the unit 
cell defined by a! and a 2 . That is, the non-commutativity 
is a result of the Aharonov-Bohm phase. T ai commutes 
with T a2 only if the flux (f> is an integer multiple of the 
flux quantum 4>q = e/h. 

If the magnetic flux <f> enclosed by a plaquette is 
(p/q)(po, where p and q are co-prime integers, then T qai 
would commute with T &2 (see Fig. 12). The simultaneous 



(8.11) 

respectively. These states are not equivalent, but have 
the same energy as i/vik since [H, T a J = 0. Therefore, 
the MBZ has a g-fold degeneracy along the b 2 direc- 
tion. Each repetition unit in the MBZ is sometimes 
called a reduced magnetic Brillouin zone. More discus- 
sions on the magnetic translation group can be found in 
Zak (1964a,b,c). 



B. Basics of magnetic Bloch band 

In this subsection, we review some basic properties of 
the magnetic Bloch band. This includes the pattern of 
band splitting due to a quantizing magnetic field, the 
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phase of the magnetic Bloch state in k-space and its con- 
nection with the Hall conductance. 

The rules of band splitting are simple in two op- 
posite limits, which are characterized by the relative 
strength between the lattice potential and the mag- 
netic field. When the lattice potential is much stronger 
than the magnetic field, it is more appropriate to start 
with the zero-field Bloch band as a reference. It was 
found that, if each plaquette encloses a magnetic flux 
(p/q)4>o, then each Bloch band would split to q sub- 
bands (Hatsugai and Kohmoto, 1990; Kohmoto, 1989). 
We know that if N is the total number of lattice sites on 
the two dimensional plane, then the number of allowed 
k-states in the Brillouin zone (and in one Bloch band) is 
N. Since the area of the MBZ (and the number of states 
within) is smaller by a factor of q, each MBB has N/q 
states, sharing the number of states of the original Bloch 
band equally. 

On the other hand, if the magnetic field is much 
stronger than the lattice potential, then one should start 
from the Landau level as a reference. In this case, if 
each plaquette has a magnetic flux (f> = (p/q)<fio, then 
after turning on the lattice potential, each LL will split 
to p subbands. The state counting is quite different from 
the previous case: The degeneracy of the original LL is 
<E>/0 n = Np/q, where <& = N<j> is the total magnetic flux 
through the two dimensional sample. Therefore, after 
splitting, each MBB again has only N/q states, the num- 
ber of states in a MBZ. 

Between the two limits, when the magnetic field is nei- 
ther very strong nor very weak, the band splitting does 
not follow a simple pattern. When the field is tuned from 
weak to strong, the subbands will split, merge, and inter- 
act with each other in a complicated manner, such that 
in the end there are only p subbands in the strong-field 
limit. 

According to Laughlin's gauge-invariance argument 
(Laughlin, 1981), each of the isolated magnetic Bloch 
band carries a quantized Hall conductivity (see Sec. II. B 
and Sec. III.C). This is closely related to the non- 
trivial topological property of the magnetic Bloch state 
(Kohmoto, 1985; Morandi, 1988). Furthermore, the dis- 
tribution of Hall conductivities among the split sub- 
bands follows a very simple rule first discovered by 
Thouless et al. (1982) et al. This rule can be derived 
with the help of the magnetic translation symmetry 
(Dana et al, 1985). We show the derivation below fol- 
lowing Dana et aZ's analysis since it reveals the important 
role played by the Berry phase in the magnetic Bloch 
state. 

In general, the phases of Bloch states at different k's 
are unrelated and can be defined independently. How- 
ever, the same does not apply to a MBZ. For one thing, 
the phase has to be non-integrable in order to account for 
the Hall conductivity. One way to assign the phase of the 
MBS itk(r) is by imposing the parallel-transport condi- 
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FIG. 13 (color online) (a) The phases of the MBS in the 
reduced MBZ can be assigned using the parallel transport 
conditions, first along the fci-axis, then along the paths par- 
alllel to the fe-axis. (b) Hyperorbits in a reduced MBZ. Their 
sizes are quantized following the Bohr-Sommerfeld quantiza- 
tion condition. The orbit enclosing the largest area is indi- 
cated by solid lines. 



tion (see Thouless's article in Prangc and Girvin (1987)), 

Ukl ° Ukl0 } = °' ( ' 8 ' 12 ' ) 

" klk2 ~&h2 Mfclfc2 ) = °' ( 8 ' 13 ) 

The first equation defines the phase of the states on the 
fci-axis; the second equation defines the phase along a line 
with fixed k\ (see Fig. 13(a)). As a result, the phases of 
any two states in the MBZ have a definite relation. 

The states on opposite sides of the MBZ boundaries 
represent the same physical state. Therefore, they can 
only differ by a phase factor. Following Eqs. (8.12) and 
(8.13), we have 



u k 1 +b 1 /q,k 2 — U kl ,k 2 '' 



.14) 
.15) 



where b\ and 62 are the lengths of the primitive vectors 
reciprocal to ai,a2. That is, the states on the opposite 
sides of the k\ boundaries have the same phase. The same 
cannot also be true for the boundaries, otherwise the 
topology will be too trivial to accommodate the quantum 
Hall conductivity. 

Periodicity of the MBZ requires that 



8{ki + bi/q) — 5(k{) + 2tt x integer. 



.16) 
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In order for the integral (l/2n) § gMBZ dk • «4.(k) (which is 
nonzero only along the upper horizontal boundary) to be 
the Hall conductivity <j h (in units of h/e 2 ), the integer 
in Eq. (8.16) obviously has to be equal to <th- 

Following the periodicity condition in Eq. (8.16), it is 
possible to assign the phase in the form, 

S(ki) = 5(h) + <T H faqai, (8.17) 

where 5(k\ + b\/q) — S(ki). On the other hand, from 
the discussion at the end of the previous subsection, we 
know that 

T ai Uk 1 k 2 = e j6, ^ 1 ^u fclfe2+ 2 7r p/ g a 2 - (8.18) 
Again from the periodicity of the MBZ, one has 

9(ki + h/q) = 6{kx) + 27rm;m e Z, (8.19) 
which gives 

0(ki) = 9(k\) + mkiqai. (8.20) 
Choose S(ki) and 9{ki) to be zero, one finally gets 

rp iqmk 1 qa 1 

- L qa 1 a k 1 k 2 — e u fci/c2+27rp/a 2 

= e i qm krqa le ip<T H k iqaiUkik ^ ( 821 ) 

But this state should also be equal to e tqkiai Uk 1 k 2 ■ There- 
fore, the Hall conductivity should satisfy 

pa H + qm=l. (8.22) 

This equation determines the Hall conductivity (mod q) 
of a MBB (Dana et ai, 1985). In Sec. VIII.D, we will see 
that the semiclassical analysis can also help us finding out 
the Hall conductivity of a MBB. 

C. Semiclassical picture: hyperorbits 

When a weak magnetic field is applied to a Bloch band, 
the electron experiences a Lorentz force and executes a 
cyclotron motion on the surface of the Fermi sea. In 
the case of the MBB, the magnetic field Bo changes the 
band structure itself. On the other hand, the magnetic 
quasi-momentum hk is a good quantum number with 
fik = 0. Therefore, there is no cyclotron motion of k 
(even though there is a magnetic field Bo ). Similar to 
the case of the Bloch band, one can construct a wave 
packet out of the magnetic Bloch states, and study its 
motion in both the r and the k space when it is subject to 
an additional weak electromagnetic field E and <5B. The 
semiclassical equations of motion that is valid under the 
one-band approximation have exactly the same form as 
Eq. (5.8). One simply needs to reinterpret ~hk, Eo(k), and 
B in Eq. (5.8) as the magnetic momentum, the magnetic 
band energy, and the extra magnetic field <5B respectively 
(Chang and Niu, 1995, 1996). As a result, when SB is not 



zero, there exists similar circulating motion in the MBB. 
This type of orbit will be called "hyperorbit" . 

Let us first consider the case without the electric field 
(the case with both E and <5B will be considered in the 
next subsection). By combining the following two equa- 
tions of motion (Cf. Eq. (5.8)), 

fik = -er x SB; (8.23) 

OP 

hr = — - hk x Q, (8.24) 
ok 

one has, 

where «(k) = 1 + £l(k)SBe/Ti. This determines the 
k-orbit moving along a path with constant E(k) = 
Eo(k) — M(k) • SB, which is the magnetic Bloch band en- 
ergy shifted by the magnetization energy. Similar to the 
Bloch band case, it is not difficult to see from Eq. (8.23) 
that the r-orbit is simply the k-orbit rotated by ir/2 and 
(linearly) scaled by the factor fi/eSB. These orbits in the 
MBB and their real-space counterparts are the hyperor- 
bits mentioned earlier (Chambers, 1965). 

The size of a real-space hyperorbit may be very large 
(if phase coherence can be maintained during the circula- 
tion) since it is proportional to the inverse of the residual 
magnetic field SB. Furthermore, since the split magnetic 
subband is narrower and flatter than the original Bloch 
band, the electron group velocity is small. As a result, 
the frequency of the hyperorbit motion can be very low. 
Nevertheless, it is possible to detect the hyperorbit us- 
ing, for example, resonant absorption of ultrasonic wave 
or the conductance oscillation in an electron focusing de- 
vice. 

Similar to the cyclotron orbit, the hyperorbit motion 
can also be quantized using the Bohr-Sommerfeld quan- 
tization rule (see Eq. (7.3)). One only needs to bear 
in mind that k is confined to the smaller MBZ and the 
magnetic field in Eq. (7.3) should be SB. After the quan- 
tization, there can only be a finite number of hyperorbits 
in the MBZ. The area of the largest hyperorbit should be 
equal to or slightly smaller (assuming SB <C B so that 
the number of hyperorbits is large) than the area of the 
MBZ (2n/a) 2 /q (see Fig. 13(b)). For such an orbit, the 
Berry phase correction r/27r in Eq. (7.3) is very close to 
the integer Hall conductivity an of the MBB. Therefore, 
it is not difficult to see that the number of hyperorbits 
should be \l/(qStf>) + (Jh\, where Scj) = SBa 2 /cf)o is the 
residual flux per plaquette. 

Because the MBZ is q-fold degenerate (see 
Sec. VIII. A), the number of energy levels produced 
by these hyperorbits are (Chang and Niu, 1995) 

D= \l/(qS^ + a H \_ (g26) 

q 

If one further takes the tunnelling between degenerate hy- 
perorbits into account (Wilkinson, 1984a), then each en- 
ergy level will be broadened into an energy band. These 
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are the magnetic energy subbands at a finer energy scale 
compared to the original MBB. 



D. Hall conductivity of hyperorbit 

According to Laughlin's argument, each of the isolated 
subband should have its own integer Hall conductivity. 
That is, as a result of band splitting, the integer Hall 
conductivity an of the parent band is split to a distri- 
bution of integers a r (there are q of them) . The sum of 
these integers should be equal to the original Hall con- 
ductivity: an = J2r a r- There is a surprisingly simple 
way to determine this distribution using the semiclassical 
formulation: one only needs to study the response of the 
hyperorbit to an electric field. 

After adding a term — eE to Eq. (8.23), one obtains, 



% ■ A E x z 



.27) 



For a closed orbit, this is just a cyclotron motion super- 
imposed with a drift along the E x direction. After 
time average, the former does not contribute to a net 
transport. Therefore the Hall current density for a filled 
magnetic band in a clean sample is, 



d 2 k 



-ep- 



E x z 



(8.28) 



where p is the number of states in the MBZ divided by the 
sample area. Therefore, the Hall conductivity is a^ losc = 
ep/SB. If the arcal electron density of a sample is p , 
then after applying a flux <fi — p/q per plaquette, the 
MBZ shrinks by q times and p = pa/q. 

How can one be sure that both the degeneracy in 
Eq. (8.26) and the Hall conductivity <t£ 1osc are integers? 
This is closely related to the following question: How 
does one divide an uniform magnetic field B into the 
quantizing part Bq and the perturbation 5B1 The proper 
way to separate them was first proposed by Azbel (1964). 
Since then, such a recipe has been used widely in the 
analysis of the Hofstadter spectrum (Hofstadter, 1976). 

One first expands the flux <f> = p/q(< 1) as a continued 
fraction, 



P 
Q 



1 



h + 



[/l: hi h, • 
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then the continued fraction is truncated to obtain various 
orders of approximate magnetic flux. For example, </>i = 

[h] = Pi/qi,h = [h,h] = P2/q2,<p3 = [fi,h,h] = 

P3/<73j •■• e tc What is special about these truncations is 
that p r /q r is the closet approximation to p/q among all 
fractions with q < q r (Khinchin, 1964). 
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FIG. 14 A parent magnetic Bloch band at magnetic field B r 
splits to D r subbands (D r = 5 here) due to a perturbation 
5B r +i . The subbands near the band edges of the parent band 
are usually originated from closed hyperorbits. The subband 
in the middle is from an open hyperorbit. 



As a reference, we show two identities that will be used 
below: 



p r +iq r 



1r 
p r q r 



fr+lQr + Qr-1, 

i-iy. 



(8.30) 
(8.31) 



According to desired accuracy, one chooses a particular 
<j) r to be the quantizing flux, and takes 8<p r = 4> r +i ~ 4>r 
as a perturbation. With the help of Eq. (8.31) one has 



QrQr+l 
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As a result, the Hall conductivity for a closed hyper- 
orbit produced by 8B r -\ = 84> r -i/a 2 is (recall that 
Pr = Pa/q r ), 



SBr-! 



= Mr 1 ? 



r-l- 
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Substitute this value back to Eq. (8.26) for D= losc (the 
number of subbands split by 6<p r ), and use Eq. (8.30), 
one has 



D 



close 



\l/(q r S<p r ) 



-close I 



r+l- 



.34) 



This is the number of subbands split from a parent band 
that is originated from a closed hyperorbit. One can see 
that the Hall conductivity and the number of splitting 
subbands are indeed integers. 

For lattices with square or triangular symmetry, there 
is one, and only one, nesting (open) hyperorbit in the 
MBZ (for example, see the diamond-shaped energy con- 
tour in Fig. 13(b)). Because of its open trajectory, the 
above analysis fails for the nesting orbit since the first 
term in Eq. (8.27) also contributes to the Hall conductiv- 
ity. However, since the total number of hyperorbits in the 
parent band can be determined by the quantization rule, 
we can easily pin down the value of <7° pcn with the help 
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of the sum rule: cr^ aicnt = ^2 r cr r . Furthermore, D° pcn 
can be calculated from Eq. (8.26) once cr° pen is known. 
Therefore, both the distribution of the ov's and the pat- 
tern of splitting can be determined entirely within the 
scmiclassical formulation. The computation in principle 
can be carried out to all orders of r. Interested read- 
ers may consult Chang and Niu (1996) and Bohm et al. 
(2003) (Chap 13) for more details. 



IX. NON-ABELIAN FORMULATION 

In previous sections, we have considered the semiclas- 
sical electron dynamics with an Abelian Berry curvature. 
Such a formalism can be extended to the cases where the 
energy bands are degenerate or nearly degenerate (eg., 
due to spin) (Culccr et al, 2005; Shindou and Imura, 
2005). Because the degenerate Bloch states have mul- 
tiple components, the Berry curvature becomes a matrix 
with non-Abelian gauge structure. We will report re- 
cent progress on re-quantizing the semiclassical theory 
that helps turning the wavepacket energy into an effec- 
tive quantum Hamiltonian (Chang and Niu, 2008). After 
citing the dynamics of the Dirac electron as an example, 
this approach is applied to semiconductor electrons with 
spin degrees of freedom. Finally, we point out that the 
effective Hamiltonian is only part of an effective theory, 
and that the variables in the effective Hamiltonian are 
often gauge-dependent and therefore cannot be physical 
varaibles. In order to have a complete effective theory, 
one also needs to identify the correct physical variables 
relevant to experiments. 



A. Non-Abelian electron wavepacket 

The wavepacket in an energy band with D-fold de- 
generacy is a superposition of the Bloch states ip nti (Cf. 
Sec. IV), 



comes, 



D 

\ w ) = zl J rf 3 i«(q.%n(q.*)lt 



(9.1) 



where ^ n |r?„(q, t)\ 2 = 1 at each q, and a(q, t) is a nor- 
malized distribution that centers at q c (i). Furthermore, 
the wavepacket is built to be localized at r c in the r- 
space. One can first obtain an effective Lagrangian for 
the wavepacket variables r c , q c , and r] n , then derive their 
dynamical equations of motion. Without going into de- 
tails, we only review primary results of such an investi- 
gation (Culcer et al., 2005). 

Similar to the non-degenerate case, there are three 
essential quantities in such a formulation. In addition 
to the Bloch energy Eo(q), there are the Berry curva- 
ture and the magnetic moment of the wavepacket (see 
Sec. IV). However, because of the spinor degree of free- 
dom, the latter two become vector-valued matrices, in- 
stead of the usual vectors. The Berry connection be- 



R-mn(q) 
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i du nc 
dq 



(9.2) 



In the rest of this section, boldfaced calligraphic fonts are 
reserved for vector-valued matrices. Therefore, the Berry 
connection in Eq. (9.2) can simply be written as H.. 
The Berry curvature is defined as, 



•^(q) = v q X It - lit X K. 



(9.3) 



Recall that the Berry connection and Berry curvature in 
the Abelian case have the same mathematical structures 
as the vector potential and the magnetic field in electro- 
magnetism. Here, 7t and !F also have the same structure 
as the gauge potential and gauge field in the non-Abelian 
SU(2) gauge theory (Wilczck and Zee, 1984). Redefining 
the spinor basis {i/) nq } amounts to a gauge transforma- 
tion. Assuming that the new basis is obtained from the 
old basis by a gauge transformation U, then TZ- and J- 
would change in the following way: 



f)Tl 

T' = XJTXJ\ 



(9.4) 



where A is the parameter of adiabatic change. 

The magnetic moment of the wavepacket can be found 
in Eq. (4.6). If the wavepacket is narrowly distributed 
around q c , then it is possible to write it as the spinor 
average of the following quantity (Culcer et al., 2005), 



M ni (q c ) 



du n \~ i dui 

(9.5) 

where H = e'*** H e iclT . That is, M = (M) = 
rfj^Ari = J2ni Vn^-niVi- Except for the extension to mul- 
tiple components, the form of the magnetic moment re- 
mains the same as its Abelian counterpart (see Eq. (4.6)). 

As a reference, we write down the equations of motion 
for the non-Abelian wavepacket (Culcer et al., 2005): 



Kk c = 
hr r . 



< F. — er c x B 
V 

,H 



(9.6) 

hk c x F, (9.7) 



ihr] = ( M ■ B - fik c ■ n) 7} 



(9.8) 



where F = (JF), and the covariant derivative 2?/2?k c = 
<9/<9k c — iR,. The semiclassical Hamiltonian inside the 
commutator in Eq. (9.7) is 

H(r c , k c ) = E Q (k c ) ~ e0(r c ) - M(k c ) • B, (9.9) 

where k c = q c + (e/7i)A(r c ). 

Equation (9.8) governs the dynamics of the spinor, 
from which we can derive the equation for the spin vector 
J, where J = (J7"), and 3 is the spin matrix, 



ifij = J7 f J,H- Kk c ■ Tl r\ 



(9.10) 
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The spin dynamics in Eq. (9.10) is influenced by the Zee- 
man energy in Ti., as it should be. However, the signifi- 
cance of the other term that is proportional to the Berry 
connection is less obvious here. Later we will see that it 
is in fact the spin-orbit coupling. 

Compared to the Abelian case in Eq. (5.8), the In- 
equation also has the electric force and the Lorentz 
force. The r c -equation is slightly more complicated: The 
derivative in the group velocity dE/dh c is replaced by 
the covariant derivative and the variables are now matri- 
ces in general. The spinor- averaged Hamiltonian matrix 
is nothing but the wavepacket energy, E = (H) . Same as 
the Abelian case, it has three terms : the Bloch energy, 
the electrostatic energy, and the magnetization energy. 
Also, the anomalous velocity in Eq. (9.7) remains essen- 
tially the same. One only needs to replace the Abelian 
Berry curvature with the non- Abelian one. 



B. Spin Hall effect 

The anomalous velocity in Eq. (9.7) that is propor- 
tional to the Berry curvature F is of great physical sig- 
nificance. We have seen earlier that it is the transverse 
current in the quantum Hall effect and the anomalous 
Hall effect (Sec. III). The latter requires spinful elec- 
trons with spin-orbit coupling and therefore the carrier 
dynamics is suitably described by Eqs. (9.6), (9.7) and 
(9.10). 

For the non-Abelian case, the Berry curvature F is 
often proportional to the spin S (see Sees. VIII. C and 
IX. E). If this is true, then in the presence of an electric 
field, the anomalous velocity is proportional to E x S. 
That is, the trajectories of spin-up and spin-down elec- 
trons are parted toward opposite directions transverse to 
the electric field. There can be a net transverse current 
if the populations of spin-up and spin-down electrons are 
different, such as the case in a ferromagnet. This then 
leads to the anomalous Hall effect. 

If the populations of different spins are equal, then the 
net electric Hall current is zero. However, the spin Hall 
current can still be nonzero. This is the source of the in- 
trinsic spin Hall effect (SHE) in semiconductors predicted 
by Murakami et al. (2003). In the original proposal, a 
four-band Luttinger model is used to describe the heavy- 
hole (HH) bands and light- hole (LH) bands. The Berry 
curvature emerges when one restricts the whole Hilbert 
space to a particular (HH or LH) subspace. As we shall 
see in Sec. IX. E, such a projection of the Hilbert space 
almost always generates a Berry curvature. Therefore, 
the SHE should be common in semiconductors or other 
materials. Indeed, intrinsic SHE has also been theoreti- 
cally predicted in metals (Guo et al. , 2008) . The analysis 
of the SHE from the semiclassical point of view can also 
be found in Culcer et al. (2005). 

In addition to the Berry curvature, impurity scatter- 
ing is another source of the (extrinsic) SHE. This is 
first predicted by Dyakonov and Pcrcl (1971a,b) (also see 



(Chazalviel, 1975)) and the same idea is later revived by 
Hirsch (1999). Because of the spin-orbit coupling be- 
tween the electron and the (spinless) impurity, the scat- 
tering amplitude is not symmetric with respect to the 
transverse direction. This is the same skew scattering 
(or Mott scattering) in AHE (see Sec. III.D.l). 

To date, most of the experimental evidences for 
the SHE belong to the extrinsic case. They 
are first observed in semiconductors (Kato et al, 
2004; Sih et al, 2005; Wunderlich et al, 2005), and 
later in metals (Kimura et al, 2007; Seki et al., 2008; 
Valcnzucla and Tinkham, 2006). The spin Hall conduc- 
tivity in metals can be detected at room temperature 
and can be several orders of magnitude larger than that 
in semiconductors. Such a large effect could be due to 
the resonant Kondo scattering from the Fe impurities 
(Guo et al., 2009). This fascinating subject is still in 
rapid progress. Complete understanding of the intrinsic 
or extrinsic SHE is crucial to future designs that would 
generate a significant amount of spin current. 



C. Quantization of electron dynamics 

In Sec. VII, we have introduced the Bohr-Sommerfeld 
quantization, which helps predicting quantized energy 
levels. Such a procedure applies to the Abelian case and 
is limited to closed orbits in phase space. In this subsec- 
tion, we report on the method of canonical quantization 
that applies to more general situations. With both the 
semiclassical theory and the method of re-quantization at 
hand, one can start from a quantum theory of general va- 
lidity (such as the Dirac theory of electrons) and descend 
to an effective quantum theory with a smaller range of 
validity. Such a procedure can be applied iteratively to 
generate a hierarchy of effective quantum theories. 

As we have mentioned in Sec. VII. D, even though a 
Hamiltonian system always admits canonical variables, 
it is not always easy to find them. In the wavepacket 
theory, the variables r c and k c have very clear physi- 
cal meaning, but they are not canonical variables. The 
canonical variables r and p accurate to linear order of 
the fields are related to the center-of-mass variables as 
follows (Chang and Niu, 2008), 

r c = r + K(n) +£(tt), 
lik c = p + eA(r)+eBx7l(7r), (9.11) 

where 7r = p + eA(r), and G a (ir) = (e/h)(lZ. x B) • 
dlZ/dTT a . This is a generalization of the Peierls substi- 
tution to the non-Abelian case. The last terms in both 
equations can be neglected in some occasions. For ex- 
ample, they will not change the force and the velocity in 
Eqs. (9.6) and (9.7). 

When expressed in the new variables, the semiclassical 
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Hamiltonian in Eq. (9.9) can be written as, 

W(r,p) = Eo(tt) -e^(r) + eE-72.(7r) 

dE - 



B 



M(tt) - ell x 



dir 



(9.12) 



where we have used the Taylor expansion and neglected 
terms nonlinear in fields. Finally, one promotes the 
canonical variables to quantum conjugate variables to 
convert 7i to an effective quantum Hamiltonian. 

Compared to the semiclassical Hamiltonian in 
Eq. (9.9), the quantum Hamiltonian has two additional 
terms from the Taylor expansion. The dipole-energy 
term eE • 1Z is originated from the shift between the 
charge center r c and the center of the canonical vari- 
able r. Although the exact form of the Berry connection 
1Z- depends on the physical model, we will show that for 
both the Dirac electron (IX. D) and the semiconductor 
electron (IX. E), the dipole term is closely related to the 
spin-orbit coupling. The correction to the Zeeman en- 
ergy is sometimes called the Yafet term, which vanishes 
near a band edge (Yafet, 1963). 

Three remarks are in order. First, the form of the 
Hamiltonian, especially the spin-orbit term and Yafet 
term, is clearly gauge dependent because of the gauge- 
dependent Berry connection. Such gauge dependence has 
prevented one from assigning a clear physical significance 
to the Yafet term. For that matter, it is also doubtful 
that the electric dipole, or the spin-orbit energy can be 
measured independently. Second, in a neighborhood of 
a fc-point, one can always choose to work within a par- 
ticular gauge. However, if the first Chern number (or 
its non-Abelian generalization) is not zero, one cannot 
choose a global gauge in which 1Z- is smooth everywhere 
in the Brillouin zone. In such a non-trivial topological 
situation, one has to work with patches of the Brillouin 
zone for a single canonical quantum theory. Third, the 
semiclassical theory based on the variables T and Ad, on 
the other hand, are gauge independent. Therefore, the 
effective quantum theory can be smooth globally. 



D. Dirac electron 

To illustrate the application of the non-Abelian 
wavepacket theory and its re-quantization, we use the 
Dirac electron as an example. The starting quantum 
Hamiltonian is 



H = cat ■ (p + eA) + f3mc 2 ■ 
= Hq + cecx ■ A — e<j)(r), 



(r) 



(9.13) 



where a and (3 are the Dirac matrices (Strange, 1998) 
and Hq is the free-particle Hamiltonian. The energy 
spectrum of Hq has positive-energy branch and negative- 
energy branch, each with two-fold degeneracy due to the 
spin. This two branches are separated by a huge energy 
gap mc 2 . One can construct a wavepacket out of the 



positive-energy eigenstates and study its dynamics un- 
der the influence of an external field. The size of the 
wavepacket can be as small as the Compton wavelength 
A c = fi/mc (but not smaller), which is two orders of mag- 
nitude smaller than the Bohr radius. Therefore, the adia- 
batic condition on the external electromagnetic field can 
be easily satisfied: the spatial variation of the potential 
only needs to be much smoother than A c . In this case, 
even the lattice potential in a solid can be considered 
as a semiclassical perturbation. Furthermore, because of 
the huge gap between branches, interbranch tunnelling 
happens (and the semiclassical theory fails) only if the 
field is so strong that electron-positron pair-production 
can no longer be ignored. 

Since the wavepacket is living on a branch with two- 
fold degeneracy, the Berry connection and curvature are 
2x2 matrices (Chang and Niu, 2008), 



^(q) 



A 



27(7 + 1) 

2 7 3 



q x cr, 



ct + X r — — q 



'7+1 



(9.14) 
(9.15) 



where 7(9) = y 1 + (hq/mc) 2 is the relativistic dilation 
factor. To calculate these quantities, we only need the 
free particle eigenstates of H (see Eqs. (9.2) and (9.3)). 
That is, the non-trivial gauge structure exists in the free 
particle already. 

It may come as a big surprise that the free wavepacket 
also possess an intrinsic magnetic moment. Straightfor- 
ward application of Eq. (9.5) gives (Chuu et al., 2009), 



M(q) 



2m7 2 



<r + K— rq 

7 + 1 



(9.16) 



This result agrees with the one calculated from the ab- 
stract spin operator S in the Dirac theory (Chuu et al., 
2009), 



Ad 



eft 



-9 



2m 7 (q) 



(W\S\W), 



(9.17) 



in which the g- factor is two. The Zeeman coupling in 
the wavepacket energy is —Ad ■ B. Therefore, this mag- 
netic moment gives the correct magnitude of the Zeeman 
energy with the correct g-factor. Recall that Eq. (9.5) 
is originated from Eq. (4.6), which is the magnetic mo- 
ment due to circulation charge current. Therefore, the 
magnetic moment here indeed is a result of the spinning 
wavepacket. 

The present approach is a revival of Uhlenbeck and 
Goudsmit's rotating sphere model of the electron spin 
but without its problem. The size of the wavepacket 
A c constructed from the positive-energy states is two or- 
ders of magnitude larger than the classical electron radius 
e 2 /mc 2 . Therefore, the wavepacket does not have to ro- 
tate faster than the speed of light to have the correct 
magnitude of spin. This semiclassical model for spin is 
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Luttinger 
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FIG. 15 (color online) Schematic plot of the semiconductor 
band structure near the fundamental gap. The wavepacket in 
the conduction band is formed from a two-component spinor. 



structure far away from this region is not essential. It is 
then a good approximation to use the k • p expansion and 
obtain the 4-band Luttinger model or the 8-band Kane 
model (Kane, 1957; Luttinger, 1951; Winkler, 2003) to 
replace the more detailed band structure (see Fig. 15). 
In this subsection, we will start from the 8-band Kane 
model and study the wavepacket dynamics in one of its 
subspace: the conduction band. It is also possible to in- 
vestigate the wavepacket dynamics in other subspaces: 
the HH-LH complex or the spin-orbit split-off band. The 
result of the latter is not reported in this review. Inter- 
ested readers can consult Chang and Niu (2008) for more 
details, including the explicit form of the Kane Hamilto- 
nian that the calculations are based upon. 

To calculate the Berry connection in Eq. (9.2), one 
needs to obtain the eigenstates of the Kane model, which 
have eight components. Similar to the positive-energy 
branch of the Dirac electron, the conduction band is 
two-fold degenerate. Detailed calculation shows that, the 
Berry connection is a 2 x 2 matrix of the form, 



certainly very pleasing since it gives a clear and heuris- 
tic picture of the electron spin. Also, one does not have 
to resort to the more complicated Foldy-Wouthuysen ap- 
proach to extract the spin from the Dirac Hamiltonian 
(Foldy and Wouthuysen, 1950). 
From the equation of motion in Eq. (9.10), one obtains 



7m 



B + E x 



Kk,,, 



(7 + l)mc 2 



(9.18) 



This is the Bargmann-Michel-Telegdi (BMT) equation 
for a relativistic electron (Bargmann et ai, 1959). More 
discussions on the equations of motion for r c and k c can 
be found in Chang and Niu (2008). 

Finally, substituting the Berry connection and the 
magnetic moment into Eq. (9.12) and using Eq(tt) — 



m 2 c , one can obtain the effective quantum 



Hamiltonian, 



H(r, p) = j(-rv)mc 2 — e</>(r) + 



Ms 



7(7 + 1) mc 2 



Mb 

7 



a B, 



x a ■ E 



(9.19) 



in which all the 7's are functions of ir and \xb — 
efr/2m. This is the relativistic Pauli Hamiltonian. At 
low velocity, 7 ~ 1, and it reduces to the more fa- 
miliar form. Notice that the spin-orbit coupling comes 
from the dipole energy term with the Berry connection, 
as we have mentioned earlier (also see Mathur (1991); 
Shankar and Mathur (1994)). 



E. Semiconductor electron 

When studying the transport properties of semicon- 
ductors, one often only focus on the carriers near the 
fundamental gap at the T-point. In this case, the band 



n 



3 



1 



1 



irxk, 



(9.20) 



where E g is the fundamental gap, A is the spin-orbit 
spit-off gap, and V = ^(S\p x \X) is a matrix element of 
the momentum operator. 

As a result, the dipole term eE • TZ becomes, 



n s 



eE • It = aE • <r x k, 



(9.21) 



where a = (eV 2 /3) [\/E 2 g - l/(E g + A) 2 ]. The coef- 
ficient a and the form of the spin-orbit coupling are 
the same as the Rashba coupling (Bychkov and Rashba, 
1984; Rashba, 1960). However, unlike the usual Rashba 
coupling that requires structural inversion asymmetry to 
generate an internal field, this term exists in a bulk semi- 
conductor with inversion symmetry but requires an ex- 
ternal field E. 

From the Berry connection, we can calculate the Berry 
curvature in Eq. (9.3) to the leading order of k as, 



2V 2 
~3~ 



1 

El 



1 



(9.22) 



In the presence of an electric field, this would generate 
the transverse velocity in Eq. (9.7), 



vt = 2eaE x (er) 



(9.23) 



As a result, spin-up and spin-down electrons move toward 
opposite directions, which results in a spin-Hall effect (see 
Sec. IX. B for related discussion). 

The wavepacket in the conduction band also sponta- 
neously rotates with respect to its own center of mass. 
To the lowest order of k, it has the magnetic moment, 



M 



eV 2 
~3h 



1 



E a + A 



(9.24) 
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With these three basic quantities, 7?., T , and M., the 
re-quantized Hamiltonian in Eq. (9.12) can be established 
as 

H(r,p) = E (Tr)-e(f)(r)+aE-aXTT+SgfiBB~, (9.25) 

where Eq includes the Zceman energy from the bare spin 
and 

AmV 2 ( 1 1 \ . „, 

In most textbooks on solid state physics, one can find 
this correction of the g-factor. However, a clear identifi- 
cation with electron's angular momentum is often lack- 
ing. In the wavepacket formulation, we see that 5g is 
indeed originated from the electron's rotating motion. 




FIG. 16 The extent of wavepacket space, parent space, and 
full space. 



F. Incompleteness of effective Hamiltonian 

Once the effective Hamiltonian Ti (r, p) is obtained, one 
can go on to study its spectra and states, without re- 
ferring back to the original Hamiltonian. Based on the 
spectra and states, any physics observables of interest can 
be calculated. These physics variables may be position, 
momentum, or other related quantities. Nevertheless, we 
would like to emphasize that, the canonical variables in 
the effective Hamiltonian may not be physical observ- 
ables. They may differ, for example, by a Berry connec- 
tion in the case of the position variable. The effective 
Hamiltonian itself is not enough for correct prediction, if 
the physical variables have not been identified properly. 

This is best illustrated using the Dirac electron as an 
example. At low velocity, the effective Pauli Hamiltonian 
is (see Eq. (9.19)), 



W(r,p) = 



— -&f> 
2m 

/i£><7 • B, 



(JL B T 



x a ■ E 



(9.27) 



which is a starting point of many solid-state calculations. 
It is considered accurate for most of the low-energy appli- 
cations in solid state. When one applies an electric field, 
then according to the Heisenberg equation of motion, the 
velocity of the electron is 



r = 



7T 

771 



Ah 



-cr x E, 



(9.28) 



where A c is the Compton wavelength. 

However, if one calculates the velocity of a Dirac elec- 
tron according to Eq. (9.7), then the result is, 



m 



2% 



x E. 



(9.29) 



That is, the transverse velocity is larger by a factor of 
two. The source of this discrepancy can be traced back 
to the difference between the two position variables: r c 



and r (see Eq. (9.11)). One should regard the equation 
for r c as the correct one since it is based on the Dirac 
theory (also see Bliokh (2005)). 

Such a discrepancy between the same physical vari- 
able in different theories can also be understood from 
the perspective of the Foldy-Wouthuysen transformation. 
The Pauli Hamiltonian can also be obtained from block- 
diagonalizing the Dirac Hamiltonian using an unitary 
transformation. Since the basis of states has been ro- 
tated, the explicit representations of all of the observ- 
ables should be changed as well. For example, the r c 
in Eq. (9.11) can be obtained by a FW rotation, fol- 
lowed by a projection to the positive-energy subspace 
(Foldy and Wouthuysen, 1950). 



G. Hierarchy structure of effective theories 

Finally, we report on the hierarchical relations for the 
basic quantities, the Berry curvature J- and the mag- 
netic moment jVt. Let us consider theories on three dif- 
ferent levels of hierarchy - I, II, and III - with progres- 
sively smaller and smaller Hilbert spaces. These spaces 
will be called the full space, the parent space, and the 
wavepacket space respectively (see Fig. 16). 

Alternative to Eqs. (9.3) and (9.5), the Berry curvature 
and the magnetic moment can be written in the following 
forms (Chang and Niu, 2008), 



(9.30) 



l£out 

ie 



2h 



^ {Eo.m — E j)H m i X R;„. (9.31) 



l£out 



where R m ; is the Berry connection, and I sums over the 
states outside of the space of interest. From Eqs. (9.30) 
and (9.31), one sees that the Berry curvature and the 
magnetic moment for theory I are zero since there is no 
state outside the full space. With the help of the states 
in the full space, one can calculate the Berry curvatures 
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and the magnetic moment in theory II and theory III. 
They are designated as (^F p ,Ai p ) and (J 7 , M) respec- 
tively. These two sets of matrices have different ranks 
since the parent space and the wavepacket space have 
different dimensions. 

On the other hand, if one starts from the parent space, 
then the Berry curvature and the magnetic moment for 
theory II is zero (instead of J- p and A4 P ). The Berry 
curvature and the magnetic moment for theory III are 
now designated as JF' and M. . They are different from 
J- and Al since the former are obtained from the sum- 
mations with more outside states from the full space. It 
is straightforward to see from Eqs. (9.30) and (9.31) that 

T = T 1 + PT P P, 

M = M' + PMpP, (9.32) 

where P is a dimension-reduction projection from the 
parent space to the wavepacket subspace. This means 
that starting from theory II, instead of theory I, as the 
parent theory, one would have the errors PT p P and 
PM. p P. On the other hand, however, whenever the 
scope of the parent theory needs to be extended, e.g, 
from II to I, instead of starting all of the calculations 
anew, one only needs additional input from JF P and Ad p 
and the accuracy can be improved easily. 

For example, in Murakami et al. (2003, 2004) 's original 
proposal of the spin Hall effect of holes, the parent space 
is the HH-LH complex. The heavy hole (or the light 
hole) acquires a non-zero Berry curvature as a result of 
the projection from this parent space to the HH band (or 
the LH band). This Berry curvature corresponds to the 
J- above. It gives rise to a spin-dependent transverse 
velocity eE x T' that is crucial to the spin Hall effect. 

Instead of the HH-LH complex, if one chooses the eight 
bands in Fig. 15 as the full space, then the Berry curva- 
tures of the heavy hole and the light hole will get new 
contributions from PT p P . The projection from the full 
space with eight bands to the HH-LH complex of four 
bands generates a Berry curvature T p = — (2V 2 /3E 2 )J~ 
(Chang and Niu, 2008), where 3 is the spin-3/2 ma- 
trix. Therefore, after further projections, we would 
get additional anomalous velocities (eV 2 /E 2 )~E x cr and 
(eV 2 /3E 2 )E x <7 for HH and LH respectively. 

X. OUTLOOK 

In most of the researches mentioned this review, the 
Berry phase and semiclassical theory are explored in the 
single-particle context. The fact that they are so useful 
and that in some of the materials the manybody effect is 
crucial naturally motivates one to extend this approach 
to manybody regime. There has also been effort to in- 
clude the Berry phase in the density functional theory 
with spin degree of freedom (Niu and Klcinman, 1998; 
Niu et al., 1999). Recently, Haldane studied the Berry 
phase and relevant quantities in the context of Fermi- 
liquid theory (Haldane, 2004). Also, the Berry curva- 



ture on the Fermi surface, if strong enough, is predicted 
to modify a repulsive interaction between electrons to 
an attractive interaction and causes pairing instability 
(Shi and Niu, 2006). In addition to the artificial mag- 
netic field generated by the monopole of Berry curva- 
ture, a slightly different Berry curvature involving the 
time-component is predicted to generate an artificial elec- 
tric field, which would affect the normalization factor 
and the transverse conductivity (Shindou and Balcnts, 
2006). This latter work has henceforth been generalized 
to multiple-band Fermi liquid with non-Abelian Berry 
phase (Shindou and Balents, 2008). Researches along 
such a path is exciting and still at its early stage. 



There has been a growing number of researches on 
the Berry phase effect in light-matter interaction. The 
Berry curvature is responsible for a transverse shift 
(side jump) of the light beam reflecting off an inter- 
face (Onoda et al, 2004a, 2006a; Sawada and Nagaosa, 
2005) . The shift is of the order of the wavelength and is 
a result of the conservation of angular momentum. The 
direction of the shift depends on the circular polariza- 
tion of the incident beam. This "optical Hall effect" 
can be seen as a rediscovery of the Imbert-Federov ef- 
fect (Federov, 1955; Imbert, 1972). More detailed study 
of the optical transport involving spin has also been car- 
ried out by Bliokh and other researchers (Bliokh, 2006a; 
Bliokh and Bliokh, 2006; Duval et al, 2006b). The sim- 
ilarity between the side jump of a light beam and analo- 
gous "jump" of an electron scattering off an impurity has 
been noticed quite early in Berger and Bergmann's re- 
view (Chien and Westgate, 1980). In fact, the side jump 
of the electromagnetic wave and the electron can be uni- 
fied using similar dynamical equations. This shows that 
the equation of motion approach being focused in this re- 
view has very general validity. Indeed, similar approach 
has also been extended to the quasiparticle dynamics in 
Bose-Einstein condensate (Zhang et al, 2006). 



Even though the Berry curvature plays a crucial role 
in the electronic structure and electron dynamics of crys- 
tals, direct measurement of such a quantity is still lack- 
ing. There does exist sporadic and indirect evidences 
of the effect of the Berry phase or the Berry curvature 
through the measurement of, for example, the quantum 
Hall conductance, the anomalous Hall effect, or the Hall 
plateau in graphene. However, this is just a beginning. 
In this review, one can see clearly that in many circum- 
stances, the Berry curvature should be as important as 
the Bloch energy. Condensed matter physicists over the 
years have compiled a huge database on the band struc- 
tures and Fermi surfaces of all kinds of materials. It is 
about time to add theoretical and experimental results of 
the Berry curvature that will deepen our understanding 
of material properties. There is still plenty of room in 
the quasi-momentum space! 
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